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ABSTRACT # 

f * 

This thesis presents a procedural sequence for processing and 
assessing the integrity of experimental strain gage data taken on piping 
in which gage rosettes are arranged in circumferential patterns at 
several sections of a piping system. The method of least squares is 
used to obtain "best” values for strain per unit loading for each gage 
element. This data is used to find principal stresses and values of 
normal and shearing stress on a normal cross section at each rosette 
location. The data is then smoothed to obtain normal and shearing stress as 
functions of angular position for each ring of rosettes. A computer 

* 

program has been devised to perform these operations and provide a 
graphical output to assess degree of fit with experimental data. Finally, 
manual computation is employed to study the equilibrium of various seg- 
ments of the piping in evaluating the overall integrity of the data. 

This procedure and program have been used in evaluating tests performed 
on the main seawater piping system of the USS Benjamin Franklin (SSBN 
640) • Some results of the dockside hydrostatic tests are presented 
herein for purposes of illustration. 
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1. Introduction, 



t> r - ORIGIN OF IDEA FOR THESIS 

The idea for this thesis was conceived while the writer worked at the 
Electric Boat Division of General Dynamics Corporation during an indus- • 
trial tour from the U. S. Naval Postgraduate School, 

While working on the analysis of strain gage data taken during dock- 
side hydrostatic tests and sea trials of the USS Casimir Pulaski (SSBN 
633 )» it became evident that a method to fully analyze the main seawater 
piping system from the strain gage data did not exist. The only tool 
available at the time was a computer program that had been developed by 
Electric Boat which determined the two principal stresses and principal 
stress axes for each individual set of rosette readings. 

The project of instrumenting the main seawater piping system on the 
USS Benjamin Franklin (SSBN 640) and taking strain gage data during dock- 
side hydrostatic tests and sea trials had been assigned to Electric Boat 
by the U. S. Navy Bureau of Ships, but no assignment had been made to 
analyze the data taken. The writer obtained the necessary permission 
from Electric Boat and the Office of the U. S. Navy Supervisor of Ship- 
building, Groton, Connecticut, to use the data in connection with a thesis 
at the U. S. Naval Postgraduate School and established liaison with 
personnel at Electric Boat connected with conducting the SSBN 64 O tests. 

Over the last nine months Mr, Carl P. Wilson, Supervisor of the 
piping section of the Marine Engineering Design Department at Electric 
Boat and his main assistant, Mr. William R. Hintz, have been most helpful 
and prompt in forwarding the test data and other related material requested. 
Without their assistance and cooperation this thesis could not have been 
possible. The writer would like to express his sincere appreciation to 
Dr. John E. Brock, Professor of Mechanical Engineering, U. S. Naval 
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Postgraduate School, for his continued patience, encouragement, and most 
capable guidance while acting as faculty advisor. 

The intention of this thesis is not only to analyze the data from the 
SSBN 640, but also to develop a procedural sequence for processing and 
assessing the integrity of strain gage data that could be utilized in the 
design and testing of future submarine seawater piping systems. 

. HOW DATA WAS OBTAINED 

Two complete sets of data were used 4 b developing this thesis. 

The first set of data was taken on 7 August, 1965 during the dockside 
hydrostatic test of the main seawater piping system of SSBN 640 in 
accordance with Electric Boat Division's Propulsion Plant Test Form 
Number 721-5 Revision A. The second set of data was taken during builder's 
sea trials on 19 and 20 September, 1965 in accordance with Electric Boat 
Division's Sea Trial Form Number 21.04 Revision A. However, for security 
reasons, no further reference is made in this theses to sea trial data. 

All strain gages utilized were BLH Type FABR-25-12S9, constantan grid, 
phenolic glass base, three element 45° superimposed grid rosettes. More 
detailed information concerning these rosettes may be obtained from BIH 
Specifications Catalog Bulletin 101 effective September 1964. 

Gages were equally spaced around a cross section of piping at 30° 
from each other. A cross section witL twelve rosettes will be referred to 
as a ring or rosette ring throughout the remainder of this thesis. The 
rosettes in each ring were oriented as shown in Figure 1,2*±« 0 Ttofsre were 
eight such rings. 

|flow 




6789 10 111212345 



Figure 1.2.1 
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The number twelve gage is located facing directly up for all rings 
except ring number 2, in which it faces port. The approximate locations 
of the rosette rings are shown in figures 1.2.2, 1.2.3 and 1.2.4 which are 
sketches of the piping system. The exact location of each ring and the 
piping system dimensions are not given. The dimensions may be found on 
Bureau of Ships Blueprints 5SBN 640 209 2117801 Rev. H and SSBN 64 O 209 
2117802 Rev. F and the exact ring locations may be obtained from either 
Electric Boat Division's Propulsion Plant Test Form Number 721-5 Rev. A 
or Sea Trial Data Form Number 21.04 Rev. A. 

Strain gage readings were taken using SR-4 Strain Indicators Type N 
manufactured by Baldwin-Lima-Hamilton Corporation. 

OBJECTIVES OF THESIS 

The basic purpose of this thesis is to develop methods of assessing 
the integrity of the data and, incidentally, to reduce the data to a form 
more convenient for subsequent analysis. To reach this goal a digital 
computer program, has been developed and employed. This program uses as 
imputs certain constants concerning the piping material and sizes and 
the raw strain gage data. It computes the two principal stresses, the 
angle measured from the direction of flow through the piping of the 
maximum principal stress, and the shear stress at each cross section and 
stress normal to each cross section. In addition, smoothed values of 
normal and shear stress every.. 01 radian around the circumference of the 
rosette rings are calculated and plots are drawn showing both smoothed and 
actual data points. Finally the program computes the three forces and 
moments at each rosette ring. 

The outputs of the computer program are analyzed to determine the 
integrity of the raw data, if the maximum principal stress exceeds yield 
strength, the agreement between smoothed and actual data values of 
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FIGURE 1.2.2 
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MAIN CONDENSER INLET PIPING 
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FIGURE 1.2.3 
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FIGURE 1.2.4 



normal and shear stress and if equilibrium between two rosette rings is 
established. 



For the following sections, the terms experimental data and 
experimental results are used to represent the actual strain gage data 
taken and the computed results from this data. 

Actual dimensions of the piping system, data, and actual numerical 
results have been eliminated from this thesis in order to keep the 
thesis unclassified. All numerical results given are the actual results 
multiplied by arbitrary constants. However, relative values of the 
results are the same as the actual results. 

. NOMENCLATURE 



a,b,c,d, *? 
«,f,g J 


Constants for smoothed value of normal stress 


E 


Young's Modulus 


e c 


Error 




Unit vectors in the ^ and ^ directions 


F(&),F,(&) 


Functions of the angle 0 


Fft* 1 

FRy > 

Fr* J 


Resultant equilibrium forces in the x,y and z directions 


F x 


Force on a cross section of the main seawater piping system 
in the direction of the ship’s bow. 


Fy 


Force on a cross section of the main seawater piping system 
in the vertical direction from ship's centerline. 


f 2 


Force on a cross section of the main seawater piping system 
in the starboard direction from ship's centerline. 


F«' 7 

FV 

F*' J 


Forces in the x 1 , j 1 , and z' directions 


i, J,k 


Unit vectors in the x, y and z directions. 


k,l,m,n, ] 


Constants for smoothed value of shearing stress 



x y7t z t x ‘ , f7 / > z/ Coordinate systems described in text 
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Mr* 


7 




Mr y 


> Resultant 


equilibrium 


Mrs 


J 




Mx 


Moment about the fore 




direction 


of the bow. 



moments in the x, y and z directions, 
and aft axis of the ship in the 



My 



Moment about the vertical axis of the ship in the upward 
direction. 



M s 



Moment about the port and starboard axis of the ship in the 
direction of starboard 



M*' 

My' 

Mb- 







Moments in the x' , y' , and z 1 directions. 



r 

t 

x 
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Radius of the piping. 

Pipe wall thickness, 
percent of test depth. 
Strain reading at any depth 
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,£2, £5 

Op 



I 



(Jrr\C!X 

OVum 



<Tn 







Matrix 

a 7 x 7 matrix. 

7x1 Matrices 

Three strains of a three element 90° rosette. 

Maximum principal stress angle 
Poisson’s ratio. 

Summation sign 
Maximum principal stress 
Minimum principal stress 

Stress normal to a cross section of piping in the direction 
of flow. 

Shear stress on the outside surface of the pipe positive in 
the direction of increasing rosette numbers. 

Shear stress in the xz plane 

Angle on piping cross section from rosette number 12. 

Stress tensor (dyadic)* 
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2, Theory^ 

THEORY FOR STRAIGHT LINE FITTING 

It is assuned that the strain gage readings will increase linearly 
with internal pressure which in turn does increase linearly with an in- 
crease in depth. This assumption is based on the following facts: 

(a) The piping material exhibits a straight line stress-strain 

curve. 

(b) The system is not loaded above the elastic limit of the 
material as verified by the maximum values of the maximum principal stress 
at every rosette. 

(c) The deflections are small as verified by the overall change 
in strain gage readings. 

(d) It is assumed, for the builder's sea trials, that hull 
deflection inputs to the piping increase linearly with depth. For dock- 
side hydrostatic tests, there was no such input. 

The other main effect acting on the system is the weight of the system. 
Since the system is always filled with seawater of nearly constant density 
the weight effect on the piping system will remain nearly constant and, 
therefore, not have any appreciable effect on the assumption of linearity. 

The method of least squares^ is used to obtain the best straight line 
connecting all data points of each element with the percentage of test 
depth taken as the independent variable and the strain gage readings 
taken as the dependent variable. The results of the method of least 
squares are given in the form of the equation: 



^Greater details of the theoretical relations are given in, or may 
easily be inferred from, corresponding subsection of section 3, which deals 
with the computer program used to perform the calculations. 

2 

* Spiegel, M.R., Schaum's Outline Series, Theory and Problems of 
Statistics (New York; Schaum. Publishing Co., 1961) p,220 
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Y = a + bx 



2 . 1.1 



where a is the zero percent test depth intercept and b is the slope of 
the best straight line connecting all data points from the surface to 
test depth. The slope, b, is also proportional to change in strain per 
percent change in test depth. The slopes are the values of the indi- 
vidual element strains, ®2> an< * ®3> use< i in the subsequent develop- 
ment of the theory. 

THEORY FOR SOLVING PRINCIPAL STRESSES AND ANGLES 
The equations defining the principal stresses in terms of the 
strains measured by a rectangular rosette are^; 

(Tm«x = ’’£.{’fiL±£aT+ _H "V(e,-e 3 y- +C2 € 2 - (eTTeTiT 1 } 2.2.1 

1-/1 1+/JI j 

CTm, N _A_V(6,-£ 3 ) l + C 2 ^-(€ (+ f 3 )] 1 ] 2.Q.-2 

2 1 i-/jl i+ A J 



where the elements of the rosette are numbered as shown in figure 2.2.1 . 

3 




Figure 2.2.1 

The direction of the maximum principal stress from element number one 
of the rosette in terms of the measured strains is: 



•0P - i_ arc +an 
r 2 



z (£±±€2} 



2 . 2.3 



3 

The derivation of equations 2,2,1, 2.2.2, and 2.2.3 is given 
in Perry and Lissner, The Strain Gage Primer , pages 136-138 
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THEORY FOR EXPERIMENTAL VALUES OF (T* 4 't* r 



Looking down on the surface of the piping, the stress picture looks 
like figure 2.3.1. 




Figure 2.3.1 

It is known that % which is perpendicular to the plane of CC and (j^ 
is equal to zero. 

To determine the stresses on a plane perpendicular to the x-y plane 
whose trace is mm, tensor theory is used. Temporarily use coordinates 
^,7^, and z to determine CTa/. } 'txy andtxg-. (The direction z is per- 
pen4icular to the plane shown in figure 2.3.1 with unit vector k .) 
From tensor theory, the stress tensor (dyadic) may be written: 



+ Gnmw Qy 6?? 2.3. J. 

<r«= T-4-I 2 32 

2.3.3 

txfr = 2.3.4- 

Since the tern J appears in equations 2.3.2, 2.3.3, and 2.3.4, solve 
for it first. 

( V$f — CTrnflK (£• 6^) <8.^ + OrviiKJ (i ) 6^ 

- OrnAx. Gds Op — (TnoiKj Sim ©p <?y 2-3-5 

Substituting equation 2.3.5 back into equations 2.3.2, 2.3.3, and 2.3.4 
results in the following equations: 



Gki = (Jmftx COS Op (<2|r * 0 " Orniio S t M Op ( * ( ) 

- Onm* COS 2 Op 4- OnniM Siw 2 -0p 2.3.6 
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2.3.7 



T*y = CTmftx. COS Op( <2^ • j ) - OYaikj SlK) -0p (e^ • J ) 

= (Term* CO S©p SliOOp - CJiniNj SlN)0p COSOp 

= /z((. Tmftx. — (jrniK) ) S/AJ 2 Op 

Tx2 = CTmftx. Cos&p (6^ • k ) — CJmiN Si A/ Op ( • k) 

= (OTmftx. COS 0~p) O “(CTrruKj SlAJ-©p) O 

= O 2.3.8 

> - THEORY FOR SMOOTHED VAUJBS OP NORMAL AMD SHEARING STRESS 

In developing the theory for the smoothed values of 0L and 'txy } 
it has been assumed that the radius of the piping has remained 
constant throughout the entire piping system. This assumption is based 
on the diameter readings taken by GD/EB of rosette rings two through 
eight. These readings indicated that the diameter of the piping around 
a rosette ring remained constant enough to use a constant value of radius 
for calculations in this thesis. 

Since a previously developed theory was not available, it seemed 
reasonable to assume that 0* and 'tKy vary only with 0 in the following 
Fourier Series; 

0n= G.+ b S 1 M <56 -A C COS 0 -h d Si SI 20 -M2 COS 20 + £ S I N 30 

+ cj cos 30 +-*•*• = F(0) 2. A * 1 

k + ISIKJ0-I- m COS 0 +■ Tl S/M 2.0 -h a COS 20 +■ p S/kj30 

-t C[ cos 30 *■•••• = F,(0) 2-4. 2 

One of the decisions that had to be made concerned the number of 
terms of the Fourier Series to use. Since curves showing the exact 
shape of the stress distribution around the circumference did not exist, 
there was no criterion available to which curves developed could be 
compared. To assist in deciding how many terms to use, the solution was 
investigated using five, seven, nine, and eleven terms. A plot of the 
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experimental values and the smoothed values of On for five, seven, nine, 
and eleven terms for rosette ring number three is shown in figure 2.4.1. 

A decision was made to use seven terms for the following reasons: 

(a) A certain amount of redundancy was necessary to permit 
good least squares evaluation. 

(b) If higher coefficients were really of significance, more 
than twelve rosettes would be required at each rosette ring to obtain 
more than seven coefficients whose accuracy could not be doubted. 

(c) For equilibrium, only the first three coefficients are used. 
Therefore, it is necessary to have these evaluated as well as possible. 

The truncated forms of equations 2.4.1 and 2.4.2 shown below will be 
used in the subsequent development of theory. 

On = a 4- bsiKJ (p + CCOS0 4- d S»W 20 +- e cos 20 f £ S/Ad 30 

- 1- <3 cos 30= F(0) 2.4.1a. 

t<Y = 1< 4 - l S I KJ 0 4 -TTl COS0 tnstw20 f O COS 20 4- pS/ AJ 30 

+ Cj COS30 = FiCtfO 2 - 4 - 2a 

The smoothed values of (fowill be discussed before proceeding with the 
smoothed values of ‘txy. 

There are N=12 stations around the circumference of the pipe for 

JB. 

which there are experimental values of On , i.e. 5 there exists a CTi 
corresponding to F(0J ), i“i,--,N. 0"i is considered to be the experi- 
mental value if F( 0 L ') is considered to be the smoothed value. It is 
assumed that: 

error = 6c » CTc — F(0c) 2.4 • 3 

It is necessary to minimize the sum of the square of the errors in order 
to have smoothed results and experimental results within reasonable 
agreement.; Therefore, 

g-Ke;) 1 
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is to be minimized. This means: 



^ = 



o 



fr° 

fro 

¥e-° 

A£ - o 

|i-0 

Equations 2.4.4 are solved in the following manner: 

ec= (Tc- F(0:) 



2.4.4 



= O'c-[a+bsiN0 + ccos#+dsiM2^4ecos20 +- fsiN 30 +<jCos30] 2.4.5 
2(e:) z » £ f<7? 2 - 2CTcF(0c) -h[F($c)] 1 } 

A£ = - 2.f C7? AHM + 2 F(0OAEiM= o 2.4. 6 

<s * § 1? 



where K = a, b, c, d, e,,f, g 

By substituting the various values of K in equation 2.4.6 the following 



equations are obtained: 

= - 2£cu + 25F(0:)=O 2.4.7 

2a 

- - 25gv sfM0: 25 F(<2>:) sinj0: = O 2.4-S 

b 

= - 22 <r; cos0<; + 25 F(0:)cos0t = O 2.4 -9 

dc 

i§.= - 25<rr siK» 20: + 22 f(0:)s/w 20 t * = o 2.4.10 

£d 

||= - 22cn:cos2<fr +■ 25 F(0O cos2< ^‘ = ° 2-4.il 
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= - 25crcsiio30c + 2S F(^)siw3f = O 2-4-ia 

5? 

=,- 2^<n cos 3</> t *+ 22 F(0 t -) cos 30r^o 2 . 4.13 



Shifting the first terms to the right hand side and substituting equation 2 - 4 - L<X. 
for F (<$l) in equations 2.4.7 through 2.4.13, the results displayed in 



matrix form are? 



M [b] - [c] 



where; 



2.4.14 




H 2s0; 5 C 

5s0i 5 s 20; fS0e0( 

5C0; ScjZjsGf; 5 C 2$ 

5s20; ls2as0- is20c0; 

5e20; £c2j&30; £c20c0; 

Ss30 c 2s3jfe0; 5s30p0,; 

5c30i 5030^0,; £c30jc0 : 



2s20; 5c20- Zs3& Sc30,- 

5s0s20; £s0p20; 2s{te30 t 5s0c30c 

Sc0s20; ic0jc20; 5 C$530; 5Lc0c30c 
Ss 2 20; 5 s20c 20; 2.s20s30; £s20p30< 

2c20p20 c 5 c 2 20 l 2c2$j 3$ Sc20c30r 

fs30s20 c 2s3JZfc20 c £ S 230 t - £s3&30c 

5c3$s20 c 5c30c20 ; Ic30s30; 5c 2 30 t 



2.4.15 



(In this matrix only, s is used for sin, c is used for cos.) 





2d; 

tG i sin0i 
l(Ticos0t 
SCFl sin20| 
2<J [ cos20f. 
S(T t sin30ii 
2O" t cos30i. 



2.4.16 



2.4.17 



20 



When the number of stations around the circumference is even and the 



stations are equally spaced, it can readily be shown that all non-diagonal 
elements of the A matrix in equation 2.4*15 become zero. In this case the 



seven equations with seven common unknowns become seven equations, each 



with only one unknown, as follows: 

a M =5c ri 

fc>S Sin l 4c = 2ov SINJ z $c 

C 2 CO5*0i - 5 gy COS l </)i 

al£sito l 2 (J); = siu* 20s \ 2.4.13 

€ 5 cos 1 2<0t = ^GVCOS V 20<‘ 

SIN 2 - 301 = 5gv S/AJ 1 30< 
cj £cos l 30J - 5 or cos 2 30 <' 

By a similar derivation, the following equations are established for 
the seven coefficients used in determining the smoothed values of Txy, 

M [bO = [c] 

where [A^ is the same as shown in equation 2.4.15* 



2.4.19 









k 

1 

m 

n 

o 

P 

q 



St* 

'It* sin0j. 
5t. ; .cos0j 
tl'i sin20r 
2't; cos20t 

{. sin30; 
It i. cos30i 



2.4.20 



2.4.21 



Equations 2.4.17 become 



kN = 

I5.sin0 = 5tisin^St 
m£coS0 = Stico £0i 
n 1 sin 2 20 ® 2tisin20t; 
o 2 cos 2 20 = $ticos20L 
pgsin^S, = Stisin 2 30‘ c 
q£ cos 2 30 — 2ticos 2 30L 



2.4.22 
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Finding the seven unknowns for On and the seven unknowns for Txyby use of 
equations 2.4.18 and 2.4.22 will be referred to as the short method in 
the remainder of this thesis. Finding the unknowns for On and ‘txyby use 
of equations 2.4.14 and 2.4.19 will be referred to as the long method in 
the remainder of this thesis. 

Once the two sets of seven unknowns have been found by either the 
short or long method, they are substituted into equations 2.4.1a and 
2.4.2a to determine the smoothed values of On and ^txy • 



THEORY FOR EQUILIBRIUM CALCULATIONS 
Equilibrium of the system was established in several steps . The 
first step was to calculate the sum of the forces and moments at each 
rosette ring cross section. The x y direction was taken as the direction of 
flow in the pipe, the y y direction as the direction from the center of 

r 

the pipe toward rosette number 12, and the z' direction as the direction 
from the center of the pipe toward rosette number 3> as shown in figure 



2.5.1. 




► x’ 



The equations to obtain the sum of the forces and moments are: 
par ptf T 

F*' = -t r On ol 0 * r"t I (a+ b siw0 +- c. cos 0 4- dsiw 20-1- • ) cl0 
=> 2iTart ° 2.5.1 



r2TT rW 

Fy= -t \ P S(U0'txy dL0 = -rtj S/w0(lc+ Isno0 + mcos0 + - • ")d0 



J o 



= - Trl r*t 



2.5.2 
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Fz ' - (j~ n r eos</> cL(/> ~ r l Coscfi (k-t /sw $ + mcostf )dL(/) 



= 77mr{~ 

'ITT 



2.5.3 



~ / C r V C ( F+/ S//U0 + MCOS</ +• “)<J$ 

Jo J Q 

= 'Zirkr^F 2.5.4- 

/Pj y ' = t r /" l S/AJ0 QajcL(^ = r v { C J/AJ0 (CL + bs/U0+CCOS$+'")cl(fo 

J Jo 



* Trir'F 

■ITT 



2.5 .s 



r^-odStfi (T/j dL<J> - rl ^J~ oo$$(<z-t F $/*></ +“)cL</> 



- - rrcrr^r 
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The second step in the equilibrium calculations is to reorient each 
individual cross section coordinate system to a standard coordinate 
system. The standard coordinate system aligns the x direction with the 
fore to aft axis of the ship, the y direction with the vertical upward 
axis of the ship^ and the z direction with the port to starboard axis of 
the ship. The reorientation is obtained using actual dimensions and 
angles taken from the system blueprints and standard coordinate re- 
orientation methods. 

The third step in the equilibrium establishment consists of making 
the segments between rosette rings into free bodies and summing the 
forces and moments acting on the segment. In addition to the forces and 
moments calculated in equations 2.5.1 through 2.5.6, the forces due to 
the fluid in the pipe segment are added to obtain the equilibrium forces 

23 



and moments. The signs of the forces and moments at the end of the 
segment in which the fluid enters are reversed because the values cal- 
culated in equations 2.5.1 through 2.5.6 are the forces of the adjacent 
segment acting in that segment and not the forces of the incoming segment 
acting in the adjacent segment as are desired for equilibrium calculations. 
Theoretically the sums of the forces and moments on each individual 
segment should be zero. 

Only step one utilizes the computer program developed. Since each 
rosette ring is oriented differently with respect to the standard 
coordinate system chosen, more time wuld be required to adjust the program 
to calculate the different reorientations for equilibrium calculations 
than to calculate the reorientations by hand. The same applies to step 
three. Therefore, steps two and three are performed by hand. 
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3. Explanation of Computer Program 

... .is ->-v . INTRODUCTORY COMMENTS 

This section presents a full explanation of the digital computer 
program developed for evaluating the strain gage data. A listing of the 
various symbols used in the program, but not defined in the program 
explanation, is found in Appendix B, A complete listing of the program 
appears in Appendix A, The statement numbers referred to in this section 
are the six digit letters following each line of the program listing.. 

The program language employed in this program is CDC FORTRAN 60. 

All computer runs were made on the Control Data Corporation 1604 computer 
and associated equipment located at the U. S. Naval Postgraduate School 
at Monterey, California. 

The program developed is divided into 12 parts. The following sub- 
sections describe these parts. 

PART 1, DATA READ IN 

Part 1 of the program reads in various quantities that remain constant 
during the entire program. These quantities most likely would vary if 
used for another piping system, but would still remain constant for that 
particular system. 

The quantities read in are N, the number of experimental depths at 
which readings were taken; E, Young's modulus for the pipe material in 
millions of pounds per square inch; ju.j Poisons Ratio for the pipe 
material used; TH(l) through TH(8), the average wall thicknesses for rings 

one through eight; RA(l) through RA(8), the pipe radius for ring one 

* 

through ring eighth; and Z, the angle between two strain gage rosettes of 
each ring in radians. Also included in this part are the commands which 

^These values of raqius and thickness were obtained by averaging 
the thicknesses and radii at the various measurement points for each ring. 
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initialize all AM, BM, and CM registers at the start of calculations 
for each ring. 

t 

Statement number 000280 indicates, the number of rings of strain gage 
rosettes for which data is to be computed. This is indicated by the 
number of times the W D0 LOOP'* is used. 

PART 2, COMPUTATION OF BEST SLOPE 

Part 2 of the program reads the actual strain gage readings in micro- 
inches per inch from data cards at the end of the program and by the 
method of least squares determines the slope of the line connecting the 
data points from 0$ to 100$ of test depth. The slope is given by the 
value of B(M) for each gage. The actual strain at zero percent depth is 
given by the calculated intercept, identified in the program by AZ. The 
zero percent depth intercept is not used in this program, however. 

An indication of the accuracy of the data can be determined in this 
part of the program. One of the quantities calculated in this part, 
identified in the program by DIFFS(I), gives the difference between each 
actual strain gage reading and the calculated best line fitting value at 
each reading depth. The magnitude of DIFFS(I) will indicate the accuracy 
of the data. 

In ’future applications DIFFS(I) could be used to expand the program 
to accept or reject each and every data point by inspecting the magnitude 
of DIFFS(I) compared to the values of DIFFS(I) for a particular gage at 
other depths. The program at present does not reject any individual data 
points. If one or more points were to be rejected, the program could 
be altered to compute a modified slope and intercept. 

The slope, B(M), in this program is given in strain per percent of 
test depth. 

Part 2 of the program is run three times before going on to part 3 of 
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the program in order to have all three strains available at once when 
computing the principal stresses for each rosette. 



PART 3, COMPUTATION OF PRINCIPAL STRESSES AND AXES 



Part three of the program utilizes, the strain per percent of test 
depth calculated in part two of the program for each element of the 
rosettes to determine the two principal stresses and direction of 
maximum principal stress taking the direction of flow within the pipe as 
the zero angle. The values of the two principal stresses are given in 
pounds per square inch per percent of test depth. 

Statements 000860, 000870, and 000880 check to see if any of the 
individual element strains are zero. If any individual strain is zero, 
the program jumps to statement 001020 which sets the two principal 
stresses and angle equal to zero. 

Statements 000890 through 001000 calculate two principal stresses 
using the formulas 



Part 2 and part 3 of the program up to statement 001050 is run twelve 
times before continuing further. The number of times these two parts are 
run is determined by the number of rosettes per rosette ring. For the 
data evaluated, there were twelve rosettes per ring. 

Statements 001060 through 001090 orient the angle of the maximum 
principal stress with respect to the direction of flow through the pipe. 
This is only necessary when one of the outside legs of the rosette is not 
oriented in the direction of flow. The rosette orientation as indicated 





and the angle using the formula 



/Wje = _L arc 2e± - (6, + e s ) 
J 2 



€,-63 
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in figure 1.2.1 has caused the necessity in this particular case for 
including statements 001060 through 001090. 

A print out of values of 0tn ftK| Gmiw and Sp for all twelve rosettes 
in a ring is called for in statements 001100 and 001110. 

PART 4 , DECISION REGARDING METHOD OF SOLUTION 
Part four of the program determines whether the short or long solution 
to solve for the smoothed values of (In and 'txy must be used. 

In subsequent sections of this thesis the use of the word inoperative 
may have one of two meanings. In one meaning, it is used to represent any 
rosette element for which data was not* recorded^. In the other meaning, 
it is used to indicate any rosette element with a zero slope as determined 
by the method of least squares. The second meaning was a decision made 
by the writer based on a close evaluation of the data without its use. 

When any element of a rosette is inoperative, the value of Onie* , 
calculated in part three, for that rosette is set equal to zero. This 
part of the program checks to see if any value of (Jm«xfor a ring is zero. 
If any one value of (Pm^x is zero, the program jumps to the long method of 
solution. If none of the values of (Tm«* for a ring is zero, the program 
goes to the short method of solution. 

PART 5, DETERMINATION OF FOURIER COEFFICIENTS - SHORT METHOD 
Part five of" the program determines the values of the coefficients 
used in the determination of the smoothed values of Cn and Txy by means 
of the short method as defined in section 2 of this thesis. The computer 
space names are used in both parts five and six in determining the same 
items. 



5When data has not been recorded for a rosette element, the 
strain readings entered on the program data card for that element are 
zero. This is done in order to keep proper data card sequence in the 
program. 
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Statement 001230 changes the angle ^>from degrees to radians while 
statements 001240 through 001260 solve for the experimental values of G"n 
and 'txy . Phi, statement number 001280, is the angle of the rosette from 
rosette number 12. Statements 001290 through 001400 compute the values 
of the sines, cosines, sines squared, and cosines squared of the single 
angle phi, twice the angle phi and three times the angle phi. The value 
of AM(l) in statement 001410 indicates the number of rosettes in each 
ring while statements 001420 through 001470 sum up the squares of the 
sines and cosines of the single, double and triple angles, phi. These 
last seven statements are the main diagonal elements of the 7x7 matrix 
shown in theory (equation 2.4.15). Statements 001480 through 001610 
calculate the values of the two 7x1 matrices shown in theory (equations 
2.4.17 & 2.4.21). Statements 001620 through 001650 change the values of 
On and ^xy into values needed in making the polar plots. The last number 
of the M D0 LOOP** which connects statements 001220 through 001660 is 
determined by the number of rosettes in the ring. 

Statements 001690 through 001780 are name changing statements. 
Statements 001790 through 001850 solve the seven equations for the 
seven unknown coefficients for the smoothed values of On while statements 
001860 through 001920 solve the seven equations for the seven unknown 
coefficients for the smoothed values of 'txy . These 14 equations are 
equations 2.4.18 and 2.4.22 in section two covering theory.. 

The remaining statements in this part are all print and format state- 
ments to give a print out of various items computed in this part. 

When part five is used in the program, statement 002010 is used to 
bypass part six of the program and go directly from part five to part 
seven. 
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PART 6, DETERMINATION OF FOURIER COEFFICIENTS - LONG METHOD 
Part six of the program determines the values of the coefficients used 
in the determination of the smoothed values of <Xn and 'txy by means of the 
long method as defined in section 2 of this thesis. 

Statement 002080 is used to determine if any rosette in the ring is 
inoperative. If any rosette is inoperative, statements 002090 through 
002150 are used and the program jpes to the end of the "DO LOOP**. This 
is used to prevent plotting an experimental value of (Tn and ‘’txy on the 
output plots when the rosette is inoperative. 

Statements 002160 through 002290 are identical to statements 001230 
through 001360 of part five. Statements 002300 through 002480 compute 
the combinations of sines and cosines of the single, double and triple 
angle phi. The items identified by AM(l) through AM(28), statements 
002490 through 002760, are the various elements of the 7x7 matrix shown 
in section two (equation 2.4.15). Statements 002770 through 002830 are 
the elements of the 7x1 matrix shown in theory (equation 2.4.17) while 
statements 002840 through 002900 are the elements of the 7x1 matrix 
shown in theory by equation 2.4.21. Statements 002910 through 002940 
are identical to statements 001620 through 001650 in part five of the 
program. The last number of the "DO LOOP" connecting statements 002070 
through 002950 is once more determined by the number of rosettes in each 
ring. 

Statements 002980 through 003460 identify all elements of the 7x7 
matrix individually (Equation 2.4.15) as can be noted by the double 
subscripts. In turn, statements 003500 through 003520 do the identical 
identification for the two 7x1 matrices (Equations 2.4.17 and 2.4.21). 
Statements 003530, 003540 and 003550 are the first three items needed in 

J} i 

the CALL GAUSS 2 statements. The two CALL GAUSS ., 2 / statements, numbers 
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003560 and 003610, call the subroutine GAUSS 2 into use. GAUSS 2 
solves for the unknowns in the matrix equations and is used twice to 
solve for the unknown coefficients of both Gw and 'uy . 

The various other statements in this part of the program are print 
and format statements which print out various quantities computed in this 
part. 

PART 7, PREPARATION FOR GRAPHICAL OUTPUT 
Part seven of the program is directly related to the subroutine DRAW, 
It is merely a series of statements which set up the labels placed on 
figures drawn by the program. More information on the titles and the 
requirement for them can be obtained by consulting Appendix D. State- 
ments 003750 through 003900 are a computed GO TO series for numbering the 
rings appropriately as computer draws the figures. 

PART 8, CONFUTATION OF SMOOTHED VALUES OF AXIAL AND SHEARING STRESS 
Part eight of the program computes the smoothed values of Gw and Tx.y 
around the entire pipe diameter and converts the values into polar 
coordinates for plotting purposes. Hie 360 degrees of the circle of the 
pipe circumference is divided into 628 segments for plotting purposes. 
Since the computer plotter used for figure drawing can only draw vertical, 
horizontal and forty-five degree lines, it was necessary to divide the 
pipe circumference into as many convenient increments as possible. The 
value 628 was chosen since the computer utilizes radians, not degrees for 
angles, and 628 is a convenient number with which to divide 27 T radians 
and 900 is the limiting number of points that can be plotted by the sub- 
routine DRAW, 

Statement 003970 divides 27T by 628. Statement 004000 computes the 
angle. Statements 004010 through 004060 compute the sines and cosines 
of the angles needed in the formulas for On and 'ho*- . Statements 004070 
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and 004080 compute the smoothed values of Gn using the formula 
< Z/=a +6s//v0 -t- ccos0 + ds/A/20 + ecos 20 y- -ps/v 30 + Cj cos 30 - 
Statements 004090 and. 004100 compute the smoothed value of '’txy using the 
formula 

txY= k-b Isi*j0 + mcos<j> + nSifJZ0 -h <0 cos 20 •/- fls/AJ 30 + jcos3<&. 

Statements 004110 through 004180 transform the values of On and txy into 
values that can be utilized in polar plotting. Statements 004110, 004120, 

i 

004150, and 004160 are used for drawing reference circles on the figures. 
These reference circles indicate the zero stress ring. Any reading inside 
the zero stress circle forO^can be taken as a compression, while any 
reading outside the zero stress circle for d N can be taken as a tension. 

For the plots of y t* l> any reading inside the zero stress circle can be 
taken as a negative shear, while any reading outside the zero stress circle 
can be taken as a positive shear. A positive shear is defined as a shear 
caused by the portion downstream of the rosette ring acting in the direction 
of increasing rosette numbers on the upstream portion. A negative shear 
is defined as a shear caused by the portion downstream of a rosette ring 
acting in the direction of decreasing rosette numbers on the upstream 
portion. The number 100, chosen for plotting the zero stress circle for 
(Jn and 40, chosen for plotting the zero stress circle for fcy, were chosen 
to keep all plotted values positive and to show as much clarity on the 
figures as possible. If, in future applications of this program, these 
scales are not convenient, change the numbers 100 and 40 in statements 
001620-001650 inclusive, 002910-002940 inclusive and 004110-004180 
inclusive to more appropriate values. 

« , . ' r • * v PART 9 , GRAPHICAL OUTPUT 

Part nine of the program consists of three statements which call the 
Subroutine DRAW for plotting the experimental and smoothed values of 
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OVj around the circumference of each pipe section at a particular rosette,^ 
a title change statement, and three more statements which call the Sub- 
routine DRAW for similar plotting of values. The CALL DRAW statements 
occur in groups of threej the first one plots a circle of zero stress, 
the second one plots the smoothed points, and the third one plots the 
experimental values. Since compression stress appears as a negative stress 
while tension stress appears as a positive stress, a polar plot could not 
be utilized unless the position of zero stress was far enough away from 
the center of the plot to prevent the stresses from being plotted 180 
degrees out of position. The proper scaling of the plots must be done 
before the program is utilized. Appendix D , the Subroutine DRAW writeup, 
gives full details on both scaling and what should be placed in a CALL 
DRAW statement. 

PART ID, CONFUTATION OF FCRCE AND MOMENT COMPONENTS 
Part ten of the program computes components of the three forces and 
the three moments at each ring which are used for the equilibrium check. 

The x-direction is taken as the direction of flow in the pipe. The y- 
direction is taken as the direction from the center of the pipe toward 
rosette number 12. The z-direction is taken as the direction from the 
center of the pipe toward rosette number three. XFOR is the force in the 
x-direction and XMOM is the moment in the x-direction. YFOR is the 
force in the y-direction and IMOM is the moment in the y-direction. 

ZFOR is the force in the z-direction and 2M0M is the moment in the z- 
direction. The three force components and three moment components are 
used in the manual computation of the three equilibrium forces and three 
equilibrium moments. 

y*, PART 11, GAUSS 2, SOLUTION OF SIMULTANEOUS EQUATIONS 

Part eleven of the program is the Subroutine GAUSS 2. This is a 
standard subroutine in the computer subroutine library at the U. S. Naval 
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Postgraduate School. The computer facility explanation of this sub- 
routine is Appendix C of this thesis, 

. PART- 12, DRAW -.GENERAL PROGRAM FOR GRAPHICAL OUTRJT 
Part twelve of the program is the Subroutine DRAW. This is also a 
standard subroutine in the computer facility subroutine library at the 
U . S. Naval Postgraduate School. The computer facility explanation of 
this subroutine is Appendix D of this thesis. 
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4. Data Evaluation 



GAGE ORIENTATION 

The gage orientation is one of the important items to which careful 
thought should be given in order to obtain the most meaningful data. The 
most important consideration concerning gage orientation is the number of 
rosettes per ring and the spacing between rosettes. As previously discussed 
in section 2, it is easier to obtain smoothed values of On and when 
an even number of rosettes equally spaced are employed because equations 
with only one unknown are more readily solved than several equations with 
several unknowns. The actual number of rosettes per ring employed in the 
SSBN 640 instrumentation appears to be a good number for the size piping 
system to be analyzed. On future instrumentation projects the size of 
the piping system should be considered in determining the even number of 
gage rosettes to be used. 

In the system analyzed the rosette orientation was shifted 45° 
alternately between even and odd numbered gages (See figure 1.2.1). 

This shift in orientation caused minor difficulties in obtaining the 
proper angle between the maximum principal stress and direction of flow 
which had to be checked very closely; it is also a source of possible 
misunderstanding. The identical results could have been obtained with all 
rosettes oriented identically and thereby eliminating these minor diffi- 
culties. It is highly recommended for future instrumentations that all 
rosettes be oriented identically. 

It might be beneficial to consider the use of T- DELTA rosettes for 
future instrumentation projects. In several rosettes, one leg indicated 
erroneous readings. If a fourth leg of the rosette had been available, 
the rosette might have produced meaningful results not otherwise obtain- 
able. 
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ACCURACY OF READINGS - DOCK TRIAL DATA 



The readings taken during dockside hydrostatic tests were recorded 
to the nearest 10y4in/in. Since it is possible to read an SR-4 Indicator, 
Type Nj more accurately than the nearest 10 ^.n/in, the accuracy of all 
readings is not considered to be the best possible. It is desirable in 
future tests where an SR-4 Indicator^Type N ; is used to read or estimate 
strain gage indications to 1 yuin/in. Seven pressure increments were 
recorded during the tests. 

One of the items printed out in Part 2 of the digital computer 
program is DIFFS(I), the difference between each actual strain gage 
reading and the least squares calculated best line fitting va]u^(j f In 
evaluating the actual data, any value of DIFFS(I) greater than 10/(Ln/in 
is considered to indicate an inaccurate reading. With this as a basis 



of inaccuracy, the following statistical tabulation of the data has been 
made. 



Table 1 



Total number of readings taken 1939 

Total number of readings with DIFFS(I) greater than 10 yuin/in 49 

$ of total in error 2.52 % 

Total number of gage elements 288 

Number of inoperative elements 11 

elements with one reading more than 10 JU' n/in in error 15 



elements with two readings more than 10 /Uin/in in error 9 
elements with three readings more than 10 yuin/in in error 4 
elements with four readings more than 10 ^Uin/in in erro r 1 



Number of 
Number of 
Number of 
Number of 

Total number of inaccurate or inoperative elements 
% of total elements that were inaccurate or inoperative 
Total number of rosettes 

Number of rosettes with at least one inaccurate reading or 
an inoperative element Ring 1 

Ring 2 
Ring 3 
Ring 4 
Ring 5 
Ring 6 
Ring 7 
Ring 8 
Total 

$ of total rosettes inaccurate or inoperative 

Total numjper of rosettes inoperative due to faulty element 
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40 
13 . 89 $ 
96 

4 

3 

4 
3 
3 
7 
2 

Jk_ 

30 

31.25$ 

9 



$ of total rosettes inoperative due to faulty element 9.375$ 

Ring 6 had one rosette with all three elements faulty 

The statistical percentages shown in Table 1 have been interpreted 

in the following manner: 

a. Faulty elements caused 9.375$ of the rosettes to be unusable 
in evaluating the system.. Although thiq percentage is considered to be 
high in obtaining good experimental results, it is not excessive in 
evaluating this system. 

b. The percentage of inaccurate rosettes caused by a reading 
error of 10yUin/in, 21. 87 5$ > does not have any real meaning because 
only one, two, or even three inaccurate readings might not have caused 
the slope of the strain depth curve to be enough in error to make the 
rosette ‘give a faulty overall stress calculation. 

c. The percentage of inaccurate readings, 2.52$, is well within 
the experimental accuracy expected, if not even less than could be expected 
from a test performed under the conditions of this test. 

PRINCIPAL STRESS RESULTS - DOCKSIDE HYDROSTATIC TESTS 
The Bureau of Ships Shipbuilding Specif ication^ MILT 16420G, states 
that the maximum allowable stress in 70-30 copper nickel piping under 
cyclic loads is 15000 psi. In publications on material properties, the 
^±eld strength of 70-30 copper nickel is 25,000 psi. The maximum value 
of (Tm**. at any rosette at 100$ test depth calculated from the dockside 
hydrostatic test data was found to be slightly less than 10,000 psi. 

This indicates that the maximum principal stress does not reach the allow- 
able stress and is well below the yield strength at any instrumented point. 
Since the rosette ring locations were chosen to be representative of the 
entire system, this is sufficient proof that in aH probability the 
stresses within the piping system do not exceed the allowable stress 
limi t of the material. It is also proof that the stress levels were 
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sufficiently low to assure elastic action, as is assumed in the analysis 
of the data. 

COMPARISON OF EXPERIMENTAL DATA TO SMOOTHED CURVES 

The sixteen figures appearing in Appendix E are polar plots of 
smoothed axial and shearing stress around the circumference of the piping 
at each rosette ring and the experimental values of the axial and shearing 
stress at the individual rosette locations. The circle ten units in 
radius on the On plots and four units in radius on Txy plots is the circle 
of zero stress on these plots. For the <Tn plots stresses inside the se w 
stress circle are Compression while those outside the circles are- 
tension*, For the plots stress inside the zero stress circle is 
shearing stress in the negative direction while stress outside the zero 
stress circle is positive shearing stress. The scales on the plots are 
arbitrary to remain within the unclassified security classification as 
noted in section 1. The stresses shown are plotted per percent of test, 
depth. 

It is readily apparent by close examination of the sixteen figures 
that the experimental values are in close agreement with the smoothed 
curves. (The experimental points are marked by an x while the smoothed 
values appear as a closed curve.) The x at the center of some plots -is 
to be ignored as it is only a necessary step in having the CDC 1604 
computer draw the plots. It will also be noticed that some plots do not 
have all twelve experimental points shown. The missing experimental 
points are an indication of the location of faulty rosettes (inoperative 
elements) as no data from fully operative rosettes was rejected in this 
evaluation. Had the missing rosettes been operative, the shape of the 
curves might have been altered. 
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EQUILIBRIUM RESULTS FROM DOCK-TRIAL DATA 
Equilibrium calculations have been performed, on the segments of piping 
between rosette rings 1 and 2, 3 and 4, and 7 and 8. Calculations were 
not performed on segments connecting any other two rings because unknown 
forces due to foundations, pipe hangers, or other constraints exist 
between these segments. Table 2 gives a tabulation on the three segments 
analyzed. 



Table 2 



SEGMENT 


Fx 


Fy 


Fz 


Mx 


My 


Mz 


Fx 


Fy 


Fz 


Mx 


My 


Mz 




Fry 


Frz Mix 


Mry 


Mrz 


1-2 


At end #1 


At end #2 


RESULTA 


NT 




-13 -53 


6 447 215 


66 


-13 


-81 


-6 


to 


116 


-65 S 


-26 j-13^To 


1259 


30 


-812 


3-4 


At end #3 


At end #4 


RESULTANT 


-10 


-14 


97 




-740 


-466 


18 


17j 


L15 


511 621 -531 


8 


jJ 


1 112|ll90 227 


3208 


7-8 


At 


* enc 


i #7 






At end 




RESULTANT 




-52 


12 


20 


-2S5 


-176 


56 


-7 


8 


-88 


1-255 


274 


122 


-59 


20 


-68 


-5 Lc | L 294 


327 



(Resultant moments are about lower numbered end of segments) 

Theoretically, for equilibrium to be satisfied, all values of Fr*., 



Fry> M(? t > and. Me? should be zero. Although only one of the 

eighteen resultant values is zero, the other values are not excessively 
high. Any one of the following influences .could have caused the non-zero 
values, 

a. The distances between rosette rings that were used in 
calculating the resultant moments were distances measured from the piping 
blueprints using the approximate locations of the rosette rings. Since 
exact values of these distances were not available, these were the best 
values available to use even though they could be a source of error in the 
calculations. 

b. The values of pipe radius and thickness used in calculating 
the forces and moments at each rosette ring were average values calculated 
by averaging the values of thickness and radius at each rosette within the 
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ring. In addition the values of radius at each rosette reported in 
Propulsion Plant Test Fora 721-5 were obtained by measuring pipe diameters 
and dividing by two, which is not a true value of the radius considering 
that the circumference of the pipe is not a true circle. Since outside 
diameter measurements for rosette ring number 1 were not recorded, the 
value of radius of ring 1 used in force and moment calculations was 
the average of the other seven rings. 

c. The -value of the forces exerted by the fluid within the pipe 
on the inside surface of the pipe was determined using average values of 
radius and wall thickness. Using average values instead of exact values 
could have once more caused errors. 

d. Errors in the Fourier coefficients caused by errors in strain 
gage readings could also have contributed to errors in the forces and 
moments calculated at each rosette ring. 

Considering all the possible sources of error in the equilibrium 
calculations, it is reasonable to expect that all the resultant forces and 
moments would not be zero. However, some of the values in table 2 are 
much larger than could be considered reasonable for the method used. 
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5 • Conclusion 



• . , - < VALIDITY CF RESULTS 

As indicated by the results shown in this thesis, the method 
developed can readily be used to obtain stress information concerning 
the system. However, the validity of the equilibrium results will depend 
largely on the accuracy of the strain gage data and other data needed for 
the equilibrium calculations. 

In order to increase the accuracy of the data and thereby increase 
the validity of the results, the following recommendations are submitted: 

a. Appoint a project engineer with experience in piping theory 
and strain gage measurements to have overall control of the instrumentation 
of the system, data taking during all tests, and reporting of results 
obtained from the data. This one person would supervise the project from 
the start of instrumentation to the submission of the final report. Having 
one man in charge of the project from start to finish will greatly reduce 
the chances of faulty communication between people working on the project, 
reduce the chances of necessary data items being omitted from the data 
taken, and reduce the chances for error in the actual data taken. 

b. Utilize automatic scanning and recording equipment to record all 

•IT* 

strain gage data during actual tests. Automatic recording equipment can 
observe strain gage reading much more accurately than equipment read 
manually and in much less time. In addition, human recording errors could 
be eliminated. A rapid record of strain gage reading at any depth could 
also be taken using the automatic scanning and recording equipment. 

c. An important aid to the person analyzing the experimental 
data from these tests would be notes taken during the period that data 
is being taken. Some of the information that would be useful as notes 
is any unusual conditions that occurred during data taking, comments 
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concerning any data that might give an immediate indication as to being 
inaccurate or unreliable, the conditions that existed during the periods 
between the taking of data,,. etc,* The project engineer should be 
observing all tests to take his own notes as well as having notes taken 
by all technicians involved in the tests, 

d. Additional data to be taken during the tests should include 
actual measured distances between two adjacent rosette rings for use in 
more accurate equilibrium calculations and the location of the welded 
seams for each section of pipe. 

By further refinement of data taking methods and actual data taken, 
the method of evaluation developed in this thesis can be used as a tool 
in analyzing any piping system instrumented in a manner similiar to the 
Main Seawater Piping System on SSBN 640. 

SUGGESTED FUTURE REFINEMENTS FOR METHOD OF EVALUATION. 

Certain refinements to the present method of evaluation have been 
considered; however, sufficient time is not .available to permit these 
ideas to be included in this thesis. The basic content of these ideas 
is listed in the following paragraphs. 

Part 2 of the program could be expanded to include a method which 
would investigate the value of DIFFS(I). This would be a two stage 
process. The zero percent depth intercept, a; slope, b; and values of 
DIFFS(l) would be found as they are in the present program. If DIFFS(I) 
were greater than a predetermined limit value, the data which produced 
that DIFFS(I) value would be rejected and new values of aero percent depth 
intercept and slope would be found. An addition of the type proposed 
would improve the accuracy of all resultant calculations. 

The present theory might be further developed to find a method to 
make some use of strain gage rosettes which have one or two inoperative 



elements. When the center element is the only one inoperative, some useful 
information may obviously be obtained from the rosette. If the orientation 
is that of the even numbered rosettes as shown in figure 1.2.1, elements 
A and C operative and B inoperative would permit evaluating the axial 
stress. If the orientation is that of the odd numbered rosettes as shown 
in figure 1.2.1, elements B and C operative and A inoperative would permit 
evaluating a shear strfess. Beyond this, if rosette failure were other 
than as indicated above, there is probably some statistical way of making 
use of the data afforded by the operative elements. However, the theory 
of such usage remains to be worked out. 

These two refinements added to the analysis method might produce more 
meaningful results from the dockside hydrostatic test data. 
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' COMPLETE PROGRAM LISTING 
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m n 



• • JO^"* 4~ , < C C H * J 

. PUOGPA-v. STRAIN 

r* : V - C NS I ON X ( 8 ) * X^AVF ( 8 ) * Y ( 8 ) *Y^AV tr ( 8 ) , I T ( 3 ) * 01 - rc ( 3 5 * : ( ? ) , AM( .‘2 ) 

i * dm ( ii ) .am i ) *m 1 1 1 1 1 ) ,?.nu : t ! ) *c.N( n • i ) .z< < ) tx* * a , : ) ,o n i { 12 

2 » *^HI 1 ( 63C1 . 51GXX ( 12 ) *TAUX V ( 12 ) . IT ! TLE* 12 ) »S IGN7( 62 0 i *TA'-N7( 63^ ) . 5 

3 M A X ( 12 ) *SMI N( 12 ) * 4NS( 12 ) * T-iETA ' 12 ) *X C * < ( 12 ) * *3<X( 12 > * * T X*< 12 ) *YT x* 

4(12) * XSNR ( ) * Y SNR ( 63 J ) * XSN7 ' 630 ) * V 5M7 ( 630 ) *X (630, * v T \ R 1 6 5 0 ) * X 

5T" <’ (6 30 ) * YT\7 (63 - ) *TH( S ) *RA( 3 ) * XFOR ( 3 ) * Yr CR( S ) * R< 3 i * xvo;.; ( £ ) , v 

fa ( 3 ) * Z ' 4 OM ( b ) 

T h- 1 »c THU S t aRT C“ PART 1. PAR~ 1 READS I N DAT*. 

READ ? v ^ * N * !: * U 
999 FORMAT < 1I4,2F15.5) 

*” H ( 1 ) = 0 • 6 2 5 
PA ( 1 ) = 6.4C 
TH(2>= 0.634 
R A ( 2 ) =6. 390 
TH(3)- 0*620 
R4( 3 >=6.409 
TH( 4 ) = 0 • 6 3 3 
RA(4)=6.374 
T H ( 5 > = 0*428 
R A ( 5 >=6.416 
T H ( 6 ) = 0 • 640 
RA (6 >=6.354 
TH ( 7 ) = 0 .423 
R A ( 7 ) = 6 . 4 1 1 
TH(8)=0.423 
RA(8>=6.410 
DO 66 K= 1 * 8 
DO 3 I =1*28 

3 AM ( I ) =0 • 

DO 4 1=1*7 
EM ( I) = 0 . 

4 CM< I ) = 0 • 

2=3.14169/6. 

C THIS IS THE END OF PART 1. 

13 DO 16 L = 1 * 1 2 

THIS IS THE START OF PART 2. PART 2 READS THE DATA AT THE END OF THE 
PROGRAM TO DETERMINE THE STRAIN PER PERCENT OF DEPTH FOR EACH 
INDIVIDUAL GAGE. 

5 J DO 68 M= 1 *3 

READ 100* ( X ( I ) *Y( I ) * I =1 *N ) 

10.- FORMAT ( 14F5.0) 

DO 996 I = 1 * N 
X S A V E ( I ) = X ( I ) 

99fa YSAVEC I ) = Y ( I ) 

Y SUM = 0. 

SUMX = C. 

SIJMY = 0. 

SUMXY = 0. 

MJMXX = n. 
y N ^ N 

ro 30 I = 1* N 
3ljM.x = SUMX +X ( I ) 

SL:'r - SUMY-t-Y(I) 

SUi’XX = SUMXX + (X( i ) *Xt I > ; 



3 o ~ r l - 
^ o ^ n z *) 
00003.0 

r) ^ r ^ 

^ n ^ r = ^ 
~ ^ /> 2 0 .** 
000070 
00005G 

r\ r\ r c ~ 

cor* ioo 
COG 110 
OQO 120 
000130 
000140 
000 160 
000160 
000170 
000180 
0 0 0 1 9 C 
000200 
000210 
000220 
000230 
000240 
000260 
000260 
000270 
000280 
000290 
000300 
000310 
000320 
000330 
000340 
000360 
000360 
000370 
000380 
000390 
000400 
000410 
000420 
000430 
000440 
000460 
000460 
000470 
000480 
000490 
000600 
000610 
000520 
000530 
000540 
C00550 
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3° SIJ M X Y = S UM X V + ( X ( I ) «■ Y C I ) ) 

D-NOM = XN* SUMXX-SUMX^S'JMX 

ANUM = SUMY * SUM XX - S'JMX * 5U V XY 

AZ = ANUM / DFNOM 

BNUM= XN*^UMXY-SMMX*StJMY 

B { M ) = BNUM / DENOM 

DO 40 J = 1*N 

YC = AZ + B ( M ) *X(J) 

DIFF = Y ( J ) - YC 

AC Y SUM = YSUM + D I FF#D I FF 
52 = YSUM / (XN - 1. ) 

S I G = S2 /DENOM 

SRA = SOR TF ( SUMX X*S I G ) 

SR 8= SQ RTF ( XN*S I G ) 

PRINT 997 * AZ* B(M)* SRA * SRB 

997 FORMAT < 2 A HC STRAIGHT LINE FITTING// 25H USING THE FORM Y = BX + A, 
1/ AH A = F 1 A • 8 / 7H B(M) = F1A.8/ 26H STANDARD DEVIATION OF A = E 1 A 
2.8/ 29H STANDARD DEVIATION OF B ( M ) = E1A.8) 

801 SUMF I T = .0 

DO 90 I = 1 * N 

Y c I T ( I ) = B ( M ) * X ( I ) + AZ 
DIFFSI I ) = Y F I T ( I )-Y ( I ) 

90 SUMF I T = SUMF I T + DIFFS(I) 

99 3 PRINT 99 » SUMF IT* { X ( I ) *YSAVE( I ) *Y( I )*YFIT( I ) * D I FFS C I ) * 1*1 tN) 

99 FORMAT (9H0SUMF IT = E 1 2 . 5/ / / 10X 1 HX 1 AX 1HY9 X 1 OHY T RANSP0SE8XAHY F I T 1 1 X 5H 
1DIFFS/// ( 5E 1 5 • 5 ) ) 

38 CONTINUE 

C THIS IS THE END OF PART 2 AND THF START OF PART 3. THIS PART 
C DETERMINES THE PRINCIPAL STRESSES AND THE ANGLE FROM THE PIPE 
C CENTERLINE OF THE MAXIMUM PRINCIPAL STRESS FOR EACH ROSETTE. 

I F ( B ( 1 ) ) AC 1 *A1C *A01 
AO 1 I F { B ( 2 ) ) A0 2*A10*A02 
A02 I F { B ( 3 ) ) A03*A10*A03 
AO 3 AA = E/2 • 

5A = B { 1 ) + B < 3 ) 

C A = B { 1 ) - B < 3 ) 

D= 1 • +U 
F = 1 . -U 

G=2 .*B { 2 ) - BA 
H=OA/F 

P=1./D*SQRTF (CA*CA+G*G ) 

SMAX<L)=AA*(H+P) 

SMIN(L)=AA*(H-P) 

ANG { L ) = 0 • 5* ( ATANFI G/CA ) + < 3.1A1 59/2. ) *< 1 .- ( CA/ABSF ( CA ) )>)*180./3.1A 
1139 

GO TO 13 
A10 SMAX I L ) =0 • 

SM I N ( L ) -0 • 

ANG { L ) = C • 

13 CONTINUE 

DO 113 L=l*ll*2 
113 THETA(L)=ANG(L)-A3. 

DO 116 L = 2 » 12 *2 
116 THETA{ L)=ANG{L) 

PRINT 3 * ( SMAX < L ) *SMIN< L ) * THETA ( L ) *L = 1 *12 ) 

5 FORMAT ( ///6X7HSMAX ( L ) 8X7HSMiN( L ) 8X8H THETA (L)///(3E13.3) ) 



OOC 560 
000370 
000580 
000590 
000600 
000610 
000620 
000630 
0006A0 
000650 
000660 
000670 
000680 
000690 
000700 
000710 
000720 
000730 
0007A0 
000750 
000760 
000770 
000780 
000790 
000800 
000810 
000820 
000830 
000 8 AO 
000850 
000860 
000870 
000880 
000890 
000900 
000910 
000920 
000930 
0009A0 
000950 
000960 
000970 
000980 
000990 
001000 
001010 
001020 
001030 
OOIOAO 
001050 
001060 
001070 
001080 
001090 
001100 
001110 
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n n n n n n n 



THIS IS T He END OF PART 3 AND Trie START OF PART A. PART 4 DETERMINES 001120 
WHETHER THE SHORT OR LONG SOLUTION TO SOLVE the STRESS DiSTRiauTIO.N 001130 

MUST 3E USED. 001140 

DO 420 L=l. 12 001150 

IF(SMAXIL) ) 42C.h30.420 001160 

420 CONTINUE 001170 

THIS IS T HE END CF PA^T 4 AND THF START OF PART 5. °ART 5 DETERMINES 00113C 
TH C VALUES OF THE COEFFICIENTS US C D IN THE THEORETICAL DETERMINATION 001190 
OF NORMAL STRESS AND SHEAR STRESS BY THE SHORT OR EVEN NUMBER OF 001200 

GAGES EQUALLY SPACED METHOD. 001210 

DO 16 L= 1 » 1 2 001220 

THETA (L )= THETA (L) *3. 14159/180. 001230 

SIGXXI L)=SMAX(L)*COSF ( THETA (L) ) *COSF I THETA! L) ) + SM IN! L )*SINF ( THETA! 001240 
1L) )*SINF(THETA(L> ) 001250 

T A'JXY ( L ) = ( SMAX ( L ) - SM 1 N ( L ) ) *0 . 5* S I NF ( 2 . *THETA ( L ) ) 001260 

XL=L 001270 

PHI(L)=XL*Z 001280 

S I =S 1 NF ( PH I I L ) ) 001290 

C I =COSF ( PH I ( L ) ) 001300 

SI2=SINF< 2.*PHI (LI) 001310 

CI2 = COSF(2.*PHI (L) ) 001320 

S13=S1NF(3.*PHI (L) ) 001330 

C13=COSF(3.«PHI <L) ) 001340 

S1SQ=S1#SI 001350 

C 1 SQ=C I *C I 001360 

S2SQ=S 12*512 001370 

C2SQ=C I 2*C 1 2 001380 

S3SQ=S13*S13 001390 

C3SQ = C I 3*C I 3 001400 

AM(1)=12. 001410 

AM(6)=AM(6)+SISQ 001420 

AM! 10) =AM< 101+C1SQ 001430 

AM ( 13) =AM ( 1 3 ) +S2SQ 001440 

AM! 15)=AM( 15J+C2SQ 001450 

AM(26)=AM(26)+S3SQ 001460 

AM ( 28 ) = AM ( 28 ) +C 3SQ 001470 

BM( 1 )=BM( 1 ) +SIGXX! L) 001480 

BM(2)=BM(2)+SI*SIGXX(L) 001490 

BM(3)=BM(3)+CI*SIGXX(L) 001500 

BM ( 4 ) =BM ( 4 ) +S 1 2*S 1 GXX ( L ) 001510 

BM(5)=BM(5)+CI2*S1GXX(L) 001520 

BM(6)=BM(6)+SI3*SIGXX(L) 001530 

BM(7)=BM( 7)+CI3*SIGXX(L) 001540 

CM! 1 )=CM( 1 l+TAUXY! L) 001550 

CM(2)=CM(2)+SI*TAUXY(L) 001560 

CM(3)=CM(3)+CI*T AUXY ( L ) 001570 

CM(4 ) = CM( 4 ) +SI 2*TAUXY( L ) 001580 

CM(5)=CM(5)+CI2*TAUXY(L) 001590 

CM(6)=CM(6)+SI3*TAUXY(L) 001600 

CM! 7)=CM< 7)+CI3*TAUXY!L) 001610 

XSXX(L) =( 100.+S1GXX! L) )*CI 001620 

YSXX(L) =( 100. +S I GXX! L) )*SI 001630 

XTXY(L) = (40. + T AUXY ( L ) ) *C I 001640 

YTXY(L)=(40.+TAUXY(L) )*SI 001650 

16 CONTINUE 001660 

PRINT 7»(SIGXX(L)*TAUXY (L) » L = 1 ,12) 001670 
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7 F QP M A T ( / / / 6 Y 3 H c I G Y y 1 iw / / / < i ■_ 1 : • - ) j 

'(1*1) - A ( 1 ) 

A|2*2)=Av,(6j 
A' 3*3) =AV( i,3) 

A ( 4 * <4 ) = A v ( 13) 

A ( 5 * 5 ) = A * ( 1*3) 

A(6#6) = A*'*(26) 

A'7*7)=A:/(28) 

DO 11 1=1*7 
I *1 )=?V( 1 ) 

: 1 CM ( 1 * 1 ) = CM( I ) 

ZX ( 1 * 1 ) - F N ( 1 * 1 ) /A( 1 , 1 ) 

Z>. (2*1 ) = 6 \ ( 2 * 1 ) / A ( 2 * 2 ) 

Z X ( 3 * } )=bN; 3*1 ) /A i 3*7) 

ZX ( 4*1 ) = bN( h* 1 ) /A( a, 4 ) 

Z X ( *? * 1 ) =CIM( 5*1 ) /A( 5* 5 > 

ZX(6*I)=BN(6*i )/A(o*6) 

ZX (7*1) -i3i\i ( 7,1 ) / A( 7, 7 ) 

X v ( 1 * 1 ) = CN ( 1 * 1 ) /A ( 1 * 1 ) 

Xr ( 2 * 1 )=( N( 2 * 1 ) /A( 2*2 ) 

XY ( 3 * 1 ) = CN ( 3* 1 ) /A ( ?* 3 ) 

XY ( A * 1 ) =CN (<**!) / A ( *+,4 ) 

XY (5*1 )=CN( 5*1) /A ( 5*5) 

XY ( 6* i ) = CN( 6* 1 ) /A ( 6*6 ) 

XY ( 7 , 1 ) =CN( 7* 1 ) /A( 7,7 ) 

PRINT 196 *ZX( 1 * 1 ) * ZX ( 2 *1 ) *ZX< 3 * 1 ) *ZX ( 4* 1 ) *ZX ( 3* i ) *Z << 6* 1 ) *ZX ( 7* 1 ) 
196 F GRMAT ( / / / 1 1 rl ZX ( 1 , i ) = E14.8//11H ZX(2*1) = Ei4.8//Iln ZX(3*1) = 

1E14.8//11H Z X ( 4 * 1 ) = E14.8//11H ZX< 5*1) = E14.8//11H ZX(6*i) - El4 
2.6//11H Z X ( 7*1) = E 1 4 • 8 / / ) 

PRINT 156 *XY ( 1 * 1 ) * ( 2 *1 ) *XY ( 3 , 1 ) ,XY ( 4* 1 ) *XY ( 5 * 1 / *XX( 6 *1 ) *XY ( 7* 1 ) 

156 FORMAT ( //.'11H XY(i,i) - E14.8//11H XY(2*i) = E14.8//11H XY(3*1) = 
1E14.6//11H X Y ( h * 1 ) ^ E14.8//11H XY ( 3, 1 ) = C14.3//11H XY(6*1) = E14 
2#b//iiH X Y ( 7 * 1 ) = E 1 4 • b / / ) 

GO 7 C 4 56 

C 1 H T S ! c . T H p OF PART 5 AND THF 0 TART OF PART 6. PART 6 DETERMINES 

THE VALUES OF THE COEFFICIENTS USED IN ^HE THEORETICAL DETERMINATION 
OF THE VALUES OF THE COEFFICIENTS DS r D IN THE THEORETICAL 
DETERMINATION OF NORMAL STRfSS AND SHFAR STRESS BY THE LONG OR 
C UNEQUALLY SPACED GAGES METHOD • 

4 3V DO 46 L = i ♦ 12 

IK(SMAXJL) ) 440*42*440 
42 ' I G X X ( L ) ~ 0 • 

T AUXY ( L ) = 0 • 

XSX X ( L ) =0 • 

YSXX ( L ) =0 • 

XTXY ( L ) =0. 

Y T X Y ( L ) =0 • 

•SO TO 46 

44'; THFTA(L)=THFTA(L) “3.14159/180. 

SIGXX ( L ) = SMAX ( L )*COSF ( THETA ( L ) } *CGSF ( T,lCTA< L ) ) ►SMIRK L ) « r 4MF ( ThETA( 
1L) )*SINr (THFTA(L) ) 

T AUXY ( L ) = ( SMAX ( L ) -"M IN < L ) ) *0«5* c - INK ( 2 • * T mLTA ( l ) i 
XL = L 

PHI ( L ) - XL^Z 
5 I =S I NF ( PH I ( L ) ) 

C I =CGS F ( PH i (L) ) 



''01690 
3016-0 
001700 
001 710 
001720 
wv i 7 3>j 
001740 
00175 0 
001760 
001770 
001780 
'“•017 £ n 
00 1900 
001810 
001820 
001830 
001840 
001650 
001860 
001870 
001880 
001990 
C 0 i 9 0 O 
001 9 1 0 
001920 
001330 
001940 
001950 
001960 
001970 
001980 
001930 
002000 
002 C 1 0 
002020 
002030 
002040 
002050 
002060 
002070 
002080 
002 090 
002100 
002110 
002120 
002130 
002 140 
002150 
002160 
002170 
002180 
002190 
002200 
002210 
002220 
002230 
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r 1 ? = rr '.'r { m 1 '(_) ) 


■'0 2 ? ^ 


7 : 3 - " ; N r ;3,f Pn[ ( L ) ) 


r 0 Z 2 6 0 


r T ? = rOSr ( 3. T (LI ) 


^ 0 2 2 7 0 


3 I 30=3 I *3 i 


0:2230 


0 I SG = C_ I *C I 


3022?: 


SC=SI*C! 


0 02 3 C C 


5152=51*5 12 


3 C 2 3 1 0 


5IC2=Si*Cl2 


002320 


CI52=CI*SI2 


0C2330 


C I C 2 = C I * C 12 


002340 


52.30 = 51 2*31 2 


002350 


5 2 r 2 = 5 1 2 * C I 2 


002360 


C. 2 5 0 = C I 2 * C I ? 


002370 


5183=81*513 


002380 


3IC3=3I*CI? 


002390 


C I S3=C 1*5!" 


002400 


CIC3=CI*CI 3 


0C2410 


32S3 = 3 I 2*S I 3 


002420 


S2C3=5I2*CI 3 


002430 


C2 5 3 = C I 2*5i3 


002440 


C2C3 = CI 2* Cl '3 


002450 


S3SQ=S I 3*SI 3 


002460 


S3C 3 = S I 3*C I 3 


002470 


C*SQ = CI 3*c I ^ 


002480 


AM( 1 ) = AN!< 1 ) +1 . 


002490 


A M { 2 )=AM(2)+SI 


002500 


A M ( 3 ) = A M ( 3 ) + C I 


002510 


AM ( 4 ) = AM ( 4 ) +S I 2 


002520 


AM(5)=AM(5)+CI2 


002530 


AM ( 6 ) = AM ( 6 ) + S I 50 


002540 


AM ( 7 ) = AM ( 7 ) +5C 


002550 


AM (8 ) = A M ( 8 ) + 5 I S 2 


002560 


AM ( 9 ) = A M ( 9 ) + S I C 2 


002 5 70 


AM ( 10)=AM( iC)+CIS0 


002580 


AM (11) = AM ( 11 )+CIS2 


002590 


A M ( 12) = A M ( 1 2 ) +C I C 2 


002600 


AM (13) = AM ( 13 ) + 52 50 


002610 


AM { 14) =AM( 1 4 ) +S2C 2 


002620 


AM ( 1 5 ) = AM { 1 5 ) +C 2 SO 


002630 


AM( 16) =AM{ 16J+SI 3 


002640 


AM ( 1 7 ) =AM ( 17) +CI3 


002650 


AM (18) = AM ( 1 8 ) + S I S 3 


002660 


AM ( 19 ) =AM ( 1 9 ) + .3 I C3 


002670 


AM(20)=AM(20)+CIS3 


002680 


AM(21)=AM(21 ) +C I C 3 


002690 


A M ( 2 2) = AM ( 22 ) + 5 2 5 3 


002700 


AM ( 2 3 ) =AM ( 2 3 ) +S2C3 


002710 


AM { 24 ) = AM ( 24 ) +C 253 


002720 


AM( 25) =AM{ 25 )+C2C3 


002730 


AM ( 26) =AM( 26 )+53 r O 


002740 


AM (27)=AM(27)+53C3 


002750 


AM ( 2 8 ) = AM ( 2 8 ) +C 3S0 


002760 


BM ( 1 )=BM( 1 ) +5 I G X X ( L) 


002770 


BM (2 )=BM( 2 ) +SI #S IGXX ( L ) 


002780 


B V (3)=BM(3)+CI*5IGXX(L) 


002790 
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7 f 



t 


) = -•'■( ^ ) + 3 I 2 * r i G'X ( L ) 


3C2600 


? *' ( 


> ) =*.v ( 3 ) + C I 2 * : I GXX ( l ) 


302610 


PV. ( 6 ) = A ) + r I V* : I r, < x { L ) 


702 9 2 '' 


Af'(7)=HV{ n ^ * 7 I C- X X ( L ) 


302830 


r ’-M 


> “ r K *' ( 1 ) + T A J X Y { L ) 


702840 


( 


: )=CM( 2 > + 3 1*7 A JXY ( L ) 


002950 


( 


?)= 0 '( 3 )+c;*^JYVf L) 


002860 


CM(4)=CM<4)+^I2* T A JX V (L) 


C02 870 


CM{ 6 )=CM( 5 ) +CI 2*TA JX* t L ) 


302880 


CM t 6 ) =C v i { 6 ) +S I 3* T A JXY ( L ) 


OC 2 393 


GU 7 )=Cf-*.( 7 ) +C I 3*T AUXY ( L ) 


002 900 


XSXX t L) =( lOC.fSIGXXt L ) ) *C I 


002 9 iO 


Y 9 X X ( L ) = ( 1 7 *° • + 9 I GX X ( L) )*SI 


032920 


XTX Y t L ) = ( 40 • + T A'JX V ( L ) ) *C I 


002930 


YTXY(L)=(40 * + T A 1 IX Y ( L ) ) * 9 I 


002940 


C^NTIN'JF 


002950 


PRINT 77* ( S IGXX( L ) * T AljX Y ( L ) ,L =1 , 12 ) 


002960 


FORMAT ( ///6X5nMGXX 10X&HT A'JXY / / / ( 2 El 5. 5 ) 


) 002970 


A ( 1 


♦ 1)=AM( 1) 


002980 


A C I 


* 2 >=a:/,< 2 ) 


002990 


A ( i 


*3) = AM( 3 ) 


003000 


a { 1 


*4 ) = AM ( 4 ) 


003010 


A ( 1 


• 5) =AM< 5 ) 


003020 


A ( 1 


• 6 ) =AM t 1 6 ) 


003030 


A ( 1 


* 7 ) = A V; ( 17) 


003040 


At 2 


• 1 ) = A M ( 2 ) 


003050 


A { 2 


* 2 ) = A v t 6 ) 


003060 


A { 2 


,3) =AM( 7 ) 


003070 


A ( 2 


,4)=AV' ( B) 


003090 


A ( 2 


*3 ) =AV( 9 ) 


003090 


A ( 2 


*6) =AM< 18 ) 


003 100 


A { 2 


9 7 ) = AM t 19) 


003110 


A ( 3 


> 1 ) = A M ( 3 ) 


003120 


A ( 3 


» 2 ) = AM ( 7 ) 


003130 


A { 3 


» 3 ) =AM t 10) 


003140 


A { 3 


9 4 ) = AM ( 11 ) 


003150 


A ( 3 


t b ) = AM t 12 ) 


003160 


A ( 3 


96 )= AM (20) 


003170 


A { 3 


9 7 ) = A V ( 2 1 ) 


003180 


A ( A 


91 ) = AM { 4 ) 


003190 


A ( 4 


9 2 ) = A M ( 8 ) 


003200 


A ( 4 


93) =AM( 11) 


003210 


A { 4 


9 4 ) * AM { 13 ) 


003220 


A ( 4 


♦ S ) = A M ( 14) 


003230 


At 4 


96 ) =AM( 22 ) 


003240 


A { 4 


,7) =AM( 23 ) 


003250 


At 6 


9 1 ) =AM ( b ) 


003260 


A ( 5 


92) =AM(9 ) 


003270 


A ( 6 


93) = AM t 12 ) 


003280 


A { 6 


,4) = AM ( 14) 


003290 


A { b 


• 5»)=AM( 15) 


003300 


At 3 


96 ) =AM( 24) 


003310 


A ( 6 


♦ 7 ) = A M ( 2 5 ) 


003320 


A t 6 


9] ) = AM ( 16) 


003330 


At 6 


9 2 ) =AM ( 18) 


003340 


A t 6 


9 3 ) = AM ( 20 ) 


003350 
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c 

c 



A ( 6 * 4 ) = A M ( 22) 

A ( 6 » 5 ) = AM ( 24 ) 

A ( 6 * 6 ) =AM ( 26 ) 

A ( 6 , 7 ) = A X ( 2 7 ) 

A ( 7 * 1 ) = A M ( * i7 ) 

A ( / » 2 ) = AM ( 1 9 ) 

A( 7*3) =AM(21 ) 

A(7,M=AM(23) 

A(7,5)=AM(25) 

A ( 7 , 6 ) = A M ( 2 7 ) 

A ( 7 * 7 ) = AM ( 28 ) 

PRINT 299 * ( ( A ( J* I) * I =1 *7) , J = 1 , 7 ) 

2 99 FORMAT I // / 2X 6 HA (1,1 ) 4 X 6 HA (1,2 ) 4 x 6 H A ( I 
1 I , 6 ) 4X6HA ( I *7 )/// ( 7E10.4 ) ) 

DO i 1 7 1=1,7 
BN (1,1 ) - B M ( I ) 

117 CM ( 1,1) =CM ( I ) 

NN = 7 
MM = 1 

EP=l«F-9 



00336C 

0C3370 

003363 

0C3390 

003400 

003410 

003420 

003430 

003440 

OC3450 

003460 

003470 

3 )4X6riA< i ,^)4X6rtA( I ,5)4X6 mA ( 00 34 80 

003490 

003500 

003510 

003520 

003530 

003540 

003550 



CALL G AU C $2 ( NN » MM , E P , A , B N » Z X , K c R ) 00 3 560 

PRINT 799*ZX(1»1)»ZX(2*1) * Z X ( 3 » 1 ) ,ZX(4,1),2X( 5 » 1 ) ,ZX(6»1)»ZX(7»1) 00 3570 

799 FORMAT ( ///11H ZX(i,l) = E14.3//11H ZX ( 2 , 1 ) = EX4.8//11H ZX(3,1) = 003580 

1E14.8//11H ZX ( 4 , 1 ) = E14.8//11H ZX(5,1) = E14.8//11H ZX(6,1) = E14 003590 

2.8//11H ZX ( 7 , 1 ) = E 1 4 • 8 / / ) 003600 

CALL GAUSS 2 ( NN , MM ,cP,A,CN,XY,KER) 003610 

PRINT 759, XY( 1,1), XT (2*1) , X Y ( 3 , 1 ) , X Y ( 4 , 1 } » X Y ( 5 * 1 ) » X Y ( 6 , 1 ) , X Y ( 7 , 1 ) 00 3620 

759 FORMAT ( ///11H XY(1,1) = E14.8//11H XY(2,1) = E14.6//11H XY ( 3 , 1 ) = 003630 

1E14.8//11H XY (4,1) = E14.8//11H XY(5,1) = E14.8//11H XY(6,1) = E14 003640 

2.8//11H X Y ( 7 , 1 > = E 1 4 • 8 / / ) 003650 

THIS IS THE END OF PART 6 AND THE START OF PART 7. PART 7 READS IN 003660 
THE TITLES THAT ARE PUT ON THE PLOTS DRAWN BY THE PROGRAM. 003670 
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DO 107 1=1*12 




003680 


107 


I T I TLE ( I >=8H 




003690 




I T I TLE ( 1 )=8H KOCH, J 


003700 




I T I TLE ( 2) =8H. 


W • , JR 


003710 




I T I TLE ( 7 ) = 8H SIG VS. 


003720 




IT I TLE ( 8 ) =8H PHI 


003730 




IT I TLE ( 9>=8HRING NO. 


003740 




GO TO (301,302 


( ,303,3^4,30 5,306,30 7,3^8) ,i< 


003750 


301 


I T I TLE ( 1 0 ) = 8H 


ONE 


003760 




GO TO 315 




003770 


302 


I T I TLE ( 1 0 ) = 8H 


TWO 


003780 




00 TO 315 




003790 


303 


I T I TL E ( 10 ) = 8H 


THREE 


003800 




GO TO 315 




003810 


304 


I T I TLE ( 1 0 ) = 8H 


FOUR 


003820 




GO TO 315 




003830 


305 


I T I T LE ( 1 0 ) = 8 H 


FIVE 


003840 




GO TO 315 




003850 


306 


I T I TLE ( 10 ) = 8H 


SIX 


003860 




GO TO 315 




003870 


307 


I T I TLE ( 1 0 ) = 8H 


SEVEN 


003880 




GO TO 315 




003890 


306 


I T I TLE ( 1 0 ) = 8H 


EIGHT 


003900 


3 1 5 


I T I TLF ( 11 ) = 8H 


3/7/66 


003910 
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LB = 4H 
LD = 8H 

C THIS IS THE END OF PART 7 AND THE START OF PART 8. PART 8 COMPUTES 
C THE THEORETICAL VALUES OF NORMAL STRESS AND SHEAR STRESS AROUND THE 
C PIPE DIAMETER AND TRANSFORFS THEM INTO VALUES THAT ARE TO BE PLOTTED. 
ZZZ=6. 28318/628. 

DO 26 1=1*628 
Z I = I 

PHI 1 ( I) =ZI*ZZZ 
SO=S I NF ( PH I 1(1)) 

CO = COSF (PHI 1 ( I ) ) 

S02=SINF(2.*PHI 1 ( I ) ) 

C02=C0SF(2.*PHI1< I ) ) 

S03=SINF< 3.*PHI 1 ( I ) ) 

C03=C0SF ( 3. *PHI 1 ( I ) ) 

S I GN7 ( I)=ZX ( 1 *1 )+ZX( 2*1 )*SO+ZX( 3*1 ) *CO+ZX ( 4 * 1 )*S02+ZX( 5* 1 ) *C02+ZX( 
16*1 )*S03+ZX< 7*1 ) *C03 

T AUN7 ( I ) =XY (1*1 )+XY( 2*1 ) *SO+XY ( 3 * 1 ) *CO+XY ( 4 * 1 ) *S02 + XY ( 5 * 1 ) *C02 + XY ( 
16*1 ) *S03+X Y (7*1 ) *C03 
XSNR ( I )=100.*CO 
YSNR ( I ) =100. *SO 
XSN7 ( I ) =( 100.+S I GN7 ( I ) )*C0 
YSN7 ( I) =( 1 00 • +S I GN7 ( I ) )*SO 
XTNR ( I ) =40 • *CO 
YTNR ( I ) =40 • *SO 
XTN7 ( I ) = ( 40 • +T AUN7 ( I ) ) *CO 
YTN7 ( I ) = ( 40 . + T AUN7 ( I ) ) *SO 
26 CONTINUE 

C THIS IS THE END OF PART 8 AND THE START OF PART 9. PART 9 CONTAINS 
C THE STATEMENTS WHICH CALL THE SUBROUTINE TO PLOT THE STRESS 
C DISTRIBUTION PATTERN. 

CALL DRAW(628*XSNR*YSNR*1*0*LD*ITITLE*50.*50.*4*4*2*2*8*8*0* LAST ) 
CALL DRAW( 62 8 * XSN7 * YSN 7 . 2 * 0 * LB * I TITLE* 50. * 50 . * 4 * 4 * 2 * 2 * 8 * 8 *0*LAST) 
CALL DRAW ( 1 2 * X SXX * YS XX . 3 * 1 *LD * I T I T LF * 50 . , 50 . * 4 * 4 * 2 * 2 ♦ 8 * 8 * 0 * L AST ) 

I T I TLE ( 7 ) =8H TAU VS. 

CALL DRAW( 628*XTNR*YTNR*1 *0*LD* ITITLE*20. *20. *4*4*2*2*8*8*0* LAST ) 
CALL DRAW(628*XTN7*YTN7*2.0*LB* I T I T LE * 2 0 . * 20 . * 4 * 4 * 2 . 2 * 8 * 8 * 0 * LAST ) 
CALL DRAW ( 1 2 * XTXY * YT XY * 3 * 1 * LD * I T I T LE * 20 • * 20 • * 4 * 4 * 2 * 2 ♦ 8 ♦ 8 * 0 * L AST ) 

C THIS IS THE END OF PART 9 AND THE START OF PART 10. PART 10 COMPUTES 
C THE THREE FORCES AND THE THREE MOMENTS AT EACH RING WHICH ARE USED 
C FOR THE EQUILIBRIUM CHECK. 

804 XFOR(K )=2.*3. 14159*ZX( 1 . 1 I *RA ( K ) *TH ( K ) 

YFOR ( K ) =-3 • 141 59*XY ( 2*1 ) *RA ( K ) *TH ( K ) 

ZFOR(K) =3. 14159*XY (3 *1 ) *RA ( K ) *TH ( K ) 

XMOM ( K ) =2.*3.14159*XY( 1*1) *RA( K)*RA(K)*TH( K) 

YMOM( K)=3#14159*ZX(2*1 )*RA(K)*RA(K)*TH(K) 

ZMOM ( K ) =-3. 14159*ZX( 3*1 )*RA(K)*RA(K)*TH(K ) 

85 CONTINUE 

PRINT 555 * ( XFOR ( K ) * YFOR ( K ) *Z FOR ( K ) *XMOM( K) ♦ YMOM ( K ) *ZMOM< K ) * K = 1 ♦ 8 ) 
555 FORMAT ( / / /2X7HXFOR ( K ) 5X7HYFOR ( K ) 5X7HZFOR ( K ) 5X7HXMOM ( K ) 5X7H YMOM ( K ) 5 
1 X7HZMOM (K)///(6E12.5) ) 

END 

C THIS IS THE END OF PART 10 AND THE START OF PART 11. PART 11 IS THE 
C SUBROUTINE GAUSS 2. 



003920 

003930 

003940 

003950 

003960 

003970 

003980 

003990 

004000 

004010 

004020 

004030 

004040 

004050 

004060 

004070 

004080 

004090 

004100 

004110 

004120 

004130 

004140 

004150 

004160 

004170 

004180 

004190 

004200 

004210 

004220 

004230 

004240 

004250 

004260 

004270 

004280 

004290 

004300 

004310 

004320 

004330 

004340 

004350 

004360 

004370 

004380 

004390 

004400 

004410 

004420 

004430 

004440 

004450 
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SUBROUTINE G AUSS 2 ( N , M * EP * A , B * X , KER ) 00 0 10 

D I MENS I ON A(11*11)*B(11*1)*X(11»1) 00020 

NPM=N+M 0Q030 

DO 2 < = 1»M 00040 

1 = N+K 00050 

D n 2 J = 1 * N 00060 

2 A(J*I) = 8 ( J 9 K) 00070 

10 DO 34 L = 1 , N 00080 

KP=0 00090 

2=0*0 00100 

DO 12 K =L ,N 00110 

I F ( Z-ABSF ( A ( K ,L ) ) ) 11 , 12*12 00120 

11 Z = ABSF ( A ( K * L ) ) 00130 

KP=K 00140 

12 CONTINUE 00150 

IF (L“KP) 13*20*20 00160 

13 DO 14 J =L , NPM 00170 

Z=A(L,J) 00180 

A ( L , J ) =A ( KP » J ) 00190 

14 A ( KP * J ) =Z 00200 

2 0 I F ( ABSF ( A ( L * L I ) -EP ) 50 t 50 * 30 00 210 

30 I F ( L-N ) 3 1 , 4 0 , 40 00220 

31 LP 1 =L+ 1 00230 

DO 34 K =LP 1 * N 00240 

IF(A(K,L)> 32,34,32 00250 

32 RATIO=A(K,L)/A(L,L) 00260 

DO 33 J = LP 1 * NPM 00270 

33 A(K,J)=A(K,J)-RATIO*A(L,J) 00280 

34 CONTINUE 00290 

40 DO 43 I = 1 * N 00300 

I I =N + 1 - I 00310 

DO 43 J=1 *M 00320 

JPN= J+N 00330 

S=0 • 0 00340 

I F ( I I -N ) 4 1 * 43 * 43 00350 

41 IIP1=II+1 00360 

DO 42 K = I IP1 * N 00370 

42 S=S+A( I I , K ) *X ( K , J ) 00380 

43 X< I I f J ) = ( A( I I , JPNJ-S ) /A ( I I * I I ) 00390 

KER= 1 00400 

RFTURN 00410 

50 KER = 2 00420 

END 00430 

C THIS IS THE END OF PART 11 AND THE START OF PART 12. PART 12 IS THE 004460 
C SUBROUTINE DRAW. 004470 

SUBROUTINE DRAW (NUMPTS* X* Y* MODCURV* ITYPE* LABEL, 00000000 

1 ITITLE, EXSCALE, YSCALE, IXUP, IYRlGHT, 00000010 

2 MODEXAX , MODE YAX , IWIDE, IHIGH, IGRID, 00000020 

3 LAST) 30 

0 40 

0 50 

A GENERAL CURVE DRAWING AND POINT PLOTTING SUBROUTINE 0000006 

PROGRAMMER J. R. WARD 00000070 

DATF FEB. 1964, RFVISED JUNE 1965 00000080 

SYSTEM FORTRAN 60 00000090 

OUTPUT LOGICAL TAPE NUMBER 8 00000100 
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MQ J £ 

I NPUT ARGUM C N T S 
1 . NUMP T S 



2. X 



3. Y 



4. MCDCURV 



5. ITYPE 



6 . LABEL 



7. I T I TLB 



8 . EXSCALE 



9. YSCALE 



ASTERISKS C H ^ \ S E K F jw FORTS •'* \ 63 0000^11 

? 0 C* 120 

— ~ p o 0 n i 3 

rs \ L> * 

N-JMopc OF POINTS TO '-< c PLOTTED. THIS MOST ALW A v $ ^ - - 13.3 

AE AT L^AST ?* AND V.'JST NOT EXCEED 30 c O< ^OINT ^ C ^ 'H 6 n 



PLOTTING* OP 9™ cqR C’J'~VE DRAWING. 

ARRAY Or X-ORDINATES. DIMENSION AT LEAST EQUAL 
TO NUMPTS A NO NOT MCRc THAN 9P n IN CALLING 
PROGRAM. 

ARRAY OF Y-0°0 1 NATES • DIMENSION AS FOR THE 
X ARRAY IN TH C CALLING PROGRAM • 

CONTROLS TH^ NUMBER O c CURVES * AND/OR SETS 0^ 
POINTS* ON FACH GRAPH. THE COD^S AR ,; 



OCQOO 170 
COO* 180 
C 0000 190 
DC00O20O 
OOOOC 2 10 
COO* 220 
0C0C0230 
00000240 
000 230 

on^OO 260 
2 C r' n ^ 0 



0 

1 

2 

3 

lONTROLS 

0 

1 

2 

3 



OMLY ONF PLOT 
FIRST PLOT ON 



ON THIS GRA D H 
MULTI -PLOT GRAPH 



POOO^RP 

y q n 



T NT^RMFOI A TP PLOT CN M'JLT I -PLOT GP APh°*°PDP 300 



LAST PLOT ON MULTI-PLOT GRAPH. 

THE TYPE Or PLOT. THE lOUES ARE 
STRAIGHT LINES JOIN SUCCESSIVE 
(STANDARD CURVE DRAWING) 

POINTS PLOTTED WITH CROSS (X) 
POINTS PLOTTED WITH PLUS (+) 
POINTS PLOTTED WITH SQUARE 



OOOOC 310 
000 320 

00000330 
PO I NT SO COCO 340 
00000330 
00000360 
00000370 
00000380 



4 POINTS PLOTTED WITH DIAMOND 

3 POINTS PLOTTED WITH TRIANGLE 

WHEN POINTS ARE S F I NG PLOT T FD !ITYPF=1 THRU. 
THE POINTS ARE NOT CONNECTED. 



00000390 
00000400 
5 ) *00000410 
00000420 



000 430 

IF A CURVE IS BEING DRAWN ( I TYPE = 0)* LABEL IS 00000440 
A 4 -CHAR ACT FR BCD CURVE IDENTIFIER WHICH WILL 6E00000450 
REPRODUCED AT THE END OF THF CURVE. FOR EX AMP L E * 0 0000 460 
LABEL = 4H ONE . IF POINTS ARE BEING PLOTTED* O0000470 
LABEL IS AN 8 -CHARACTER IDENTIFIER. THE FIRST 4 00000480 
CHARACTERS ARE REPRODUCED WITH THE FIRST PLOT T ED00000490 
POINT* AND THE LAST 4 WITH THE LAST POINT. SET 00000500 
TO BLANK ANY UNWANTED CHARACTERS. 00000310 

000- 320 

AN ARRAY OF TWELVE 8 -CHARACTER BCD WORDS. 00000530 

THE FIRST SIX WORDS WILL BE REPRODUCED AS THE OQQ00540 
FIRST LINE OF GRAPH TITLE* AND THE LAST SIX 00000550 

WORDS WILL FORM THF SECOND LINE. TH P TITLE MUST 00000560 
INCLUDE THE USERS JOB IDENTIFICATION. DIMENSION 000005/0 
12 IN CALLING PROGRAM* AND SET TO BLANK ALL 00000580 

UNWANTED CHARACTERS. 00000590 

000 600 

X-SCALE (UNITS/INCH) AS POSITIVE FLOATING POINT 00000610 
VARIABLE WITH ONE FIGURE SIGNIFICANCE. SET TO 00000620 

ZERO FOR AUTO-SCALE. 0000C630 

000 640 

Y-SCALE (UNITS/INCH) AS POSITIVE FLOATING POINT 00000650 
VARIABLE WITH ONE EIGURF SIGNIFICANCE. SET TO 00000660 
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c 
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ZERO FOR AUTO-SCALE. 



10 . ixup 



11. IYRIGHT 



12. MODE XAX 



13. MODE Y AX 



14. IWIDF 



15. IHIGH 



16. IGRID 



17. LAST 



X-AXIS OFFSET FROM BOTTOM OF GRAPH I iM INCHES. 
THIS MUST NOT EXCEED I H I Gh * AND MUST NOT 8E 
NEGATIVE. 



graph IN 
MUST NOT 



0 0 0 0 C 6 7 0 
COO 680 
00000690 
00000700 
00000710 
000* 720 
INCHES. TH I S00000 730 
BE NEGATIVE. 00000740 
000 750 

MODE OF X-AXIS OFFSET. SEE CODES BELOW. 00000760 

000 770 

Y-AXIS OFFSET. THE CODES ARE AS FOLLOWS 00000780 
COMPUTED OFFSET* HOLDING ORIGIN ON 00000790 
GRAPH. THE CORRESPONDING IXUP OR 00000300 

IYRIGHT IS IGNORED 00000810 

COMPUTED OFFSET* WITH ORIGIN OFF THE 00000820 
GRAPH IF APPROPRIATE. THE CORR ESPOND-OOOOO 8 30 
ING IXUP OR IYRIGHT IS IGNORED. USE 00000840 
ONLY WITH AUTO-SCALE 

axis offset as specified by ixup or 

IYRIGHT. 



Y-AXIS OFFSET FROM LEFT OF 
MUST NOT EXCFED IWIDE* AND 



MODE OF 
0 



WIDTH OF GRAPH IN INCHES. THIS MUST NOT EXCEED 
NINE. ZERO WILL BE READ AS EIGHT INCHES. 



HEIGHT OF GRAPH IN INCHES. 
c I FTEEN. ZERO WILL BE READ 



THIS MUST NOT EXCEED 
AS EIGHT INCHES. 



IF SET TO 1* A ONE INCH BY ONE 
BE SUPERIMPOSED ON THE GRAPH. 



INCH GRID WILL 



INDICATES TO CALLING PROGRAM WHETHFR LAST PLOT 
WAS COMPLETED SUCCESSFULLY. THE CODES ARE 



00000850 
00000860 
00000870 
000* 880 
00000890 
00000900 
000 910 
00000920 
00000930 
000 940 
00000950 
00000960 
0000 970 
00000980 
00000990 
00001O00 
00001010 
00001020 



0 LAST PLOT COMPLETED SUCCESSFULLY 

1 LAST PLOT NOT SUCCESSFUL 

2 LAST PLOT NOT SUCCESSFUL AND NO 
FURTHER GRAPH OUTPUT WILL BE ATTEMPT-00001030 
ED UNTIL MODCURV IS NEXT ONE OR ZERO 00001040 

3 DRAW WAS ENTERED WITH MODCURV NOT 
EQUAL TO ONE OR ZERO WHILE THE ERROR 
LOCK-OUT WAS SET. 

THIS ARGUMENT MUST ALWAYS BE A NAME IN THE CALL 
STATEMENT. NEVER A NUMBER. 



NOTE 

ALL ARGUMENTS FROM NUMBER 7 THRU. NUMBER 16 ARE IGNORED WHEN 
MODCURV IS EITHER 2 OR 3. HOWEVER* ARGUMENTS MUST NEVER BE 
OMITTED FROM THE CALLING STATEMENT. IT IS MERELY THEIR VALUES 
WHICH ARE THEN IRRELEVANT. ARGUMENTS MAY BE LISTED BY NAME OR 
VALUE IN THE CALL STATEMENT. NO VALUE IN THE CALLING PROGRAM 
WILL BE ALTERED BY THIS SUBROUTINE. 

RFFFRENCF 

THE BINARY TAPE FORMAT REQUIRED BY THE OFF-LINE PLOTTER IS 
DFSCP I BED IN THE WRITFUP OF THE CDC 160A GRAPH PLOT PROGRAM 



00001050 

00001060 

00001070 

00001080 

00001090 

00001100 

00001110 

00001120 

00001130 

00001140 

00001150 

00001160 

00001170 

00001180 

00001190 

00001200 

00001210 

00001220 
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n n n n 



1 

2 



* 



9079 

9071 



1001 

1002 

1003 

1004 



1003 

1006 

1007 

1100 



1008 

1101 



1009 

1010 
1011 

100^ 

1 

1 0 0 
2 

9000 
9100 

9001 



(Tner^T, c = <F ' r ]i # T^r F 9 R A T o'- ^ ■ ; 'jcf) ay ' r uc rr,r } £ ~ v ? c 

SIMILAR c X 0 1 - T T-aT Th L i\Tr9^CLATI0\ AV'3; V EM vsi r L Z i** • 



DIMENSION X ( 9 0 0 ) , Y(VOC)* ITIt L £(12!)* J X T I T ( 1 2 ) » J Y T I T ( 1 2 J * 

LTITLE(lu), ,< A X I S ( 3 1 , I CUR V ( 4f>0 ) * JG'!l><25)* ICvNTC), 
JJTITLE( : 2 ) 

I 90 J NT = I TYP C 

CON ( ICONTRL = 4900"“ ) . 

CON ( I CiJRy 3 - 77 ‘ 7 292 rk 2 "20a , ic'JR\/4 = r j * ) # 

9 ENLACE ’a! t h D a ^ s. ST A T ? V E M r IN FORTRAN 6 2-3. 

PUT - r j ^ rj AT A STATcvcNT, 

IF (JTEST - 73346912) *0 “>0 , 9r 7 1 * 9 r 7 r 

ITEST = o 



2 U"' 

" ir “?126C 

o::oi267 

''0001270 
^ 0001299 
00 OC 1 2 q C 
0 C 0 1 3 00 
n 0 001310 
^ 9 ^ ^ 1 32 ^ 
00001330 

000^1340 

a n * r 1 330 

0000 1360 



JTEST = 73546912 
CONTINUE 

REMOVE ABOVE NONSENSE IN FORTRAN 63. 

CHECK PREVIOUS OPERATION OF ROUTINE* IF ANT. CODES ARE 
ITEST = C IF PREVIOUS GRAPH, IF ANY, COMPLETED 

ITEST =1 IF PREVIOUS GRAPH NOT COMPLETED 

ITEST = 2 IF ERROR FOUND WHILE MODCURV WAS ONE, OP IF 
M ODCURV WAS ILLEGAL. 



00091370 
o 000 1 360 
90001390 
O00014O0 
0000141Q 
00001420 
00001430 
O0C9 1440 



IFIITFST - 2 ) 1O0O, inni , 1900 
I F ( MODCURV) l n 03,l 0^2, 1003 
ITEST = 0 
GO TO 1000 
IF (MODCURV 
LAST = 3 
RETURN 

SET UP 



- 1 ) 1004, 1 002* 1004 



ERROR RETURN ROUTINE 
IF ( I TEST ) 1009, 1006 ,1009 
I F ( MODCURV) 1007 ,1003 ,1007 
PRINT HOC 

FORMAT ( 59H NO FURTHER GRAPH OUTPUT 
ONE. ,/) 

ITEST = 2 
LAST = 2 



ENTRY AT STATEMENT 1005. 



UNTIL MO DC JR V NFXT IS ZERO 



09001430 
00901469 
00001470 
00001480 
90001490 
00001 500 
00001510 
00001520 
00001530 
00001540 
00001550 
OROOOO 1 5 60 
000C1 570 
00001580 
00001590 



return 

PRINT 1101 

FORMAT ( 30H THIS PLOT WILL NOT 3E OUTPUT. ,/) 

LAST = 1 

RETURN 

IF (MODCURV - 2) 101°, 1008*1010 
IF (MODCURV - 3)1007, 1011 , 1 C C 7 
ITEST = 0 
GO TO 1098 

CHECK LEGALITY O c INPUT ARGUMENTS. 

I F ( NUMPTS -2)1 *2,2 
PRINT 100 

FORMAT (/, 32H NUMPTS MUST NOT BF LESS THAN 2. 
GO TO 1095 

IF ( I POINT ) 9 nr o,9 n n 4, 99 01 
PRINT 9100 

FORMAT (/, 15H ILLEGAL ITYPE. ) 

GO TO 1 OC 5 

IFCIPOINT - 5 ) 90n2 ,9 n 02 *9000 



00001690 
00001610 
90001620 
00001630 
00001640 
00001650 
00001660 
00001670 
00001680 
00001690 
00001700 
00001710 
00001720 
00001730 
00001740 
90001 750 
00001760 
00001770 
00001780 
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900? I F ( NU V PTS ~ 3 ri )3*3»9' s 03 

qoh- 5 DDtf.tr 9101 

910 1 FORMAT (/, 46H NU^PT r « MUST NOT SXCcp-> F q^ oq jvt PLOT T ! NO # ) 

GO TO 10 r S 

900^ IFfNUMDTS - , q , opnc 

9005 P 3 I N T 9102 

9102 FORMAT (/, 2SH NUMPTS MOST NOT EXCEED 9^0. ) 

GO TO 1005 
3 IX = 1HX 
IY = 1HY 

AMAXX = -0.2E+100 
AMAXY = -0* 2E+ 1 00 
AMINX = +0.2E+10C 
AM I NY = +O.2E+100 

1= 1 , numpts 

M AX 1 F ( X ( I ) , AMAXX) 

MAX IF ( Y ( I ) * AMAXY ) 

MINIFIX ( I ) tAMlNX) 

MIN1FIYI I ) * AM I NY) 

MAX1FI A8SFI AMAXX) » AdSF(AMAXY), ABSF(AMINX), ASSF(AMINY)) 
1.0E+99)1022tl022»1021 



102C 



58H NO X OR Y VALUF MAY EXCEED 1.0F+99 IN 



1 ,08-97) 1023, 1025,1 025 
1 .0E-O7) 1024,1025,1025 



DO 1020 
AMAXX = 

AMAXY = 

AM I MX = 

AMI NY = 

AMAXA = 

I F ( AM AX A - 

1021 PRINT 1102 

1102 FORMAT (/, 

1ITUDE. ) 

GO TO 1005 

1022 I F ( ABSFIAMAXX - AMINX) 

1023 IFIABSFI AMAXY - AMINY) 

1024 PRINT 1103 

1103 FORMAT (/, 38H ALL POINTS HAVF THE SAME COORDINATES, 
GO TO 1005 

1025 I F { I TEST ) 4 , 7 ,4 

4 I F { MODCURV - 2)5,240,5 

5 I F { MODCURV - 3)6,240,6 

6 PRINT 101 

101 FORMAT (/, 17H ILLEGAL MODCURV. ) 

GO TO 1005 

7 I c (MODCURV) 6,9,3 

8 I F( MODCURV - 1 )6,9,6 

9 I F ( IWIDF) 10,11 ,] 2 
in 1 TIT = 5HIWI0F 

PRINT 102, ITIT, ITIT 



102 FORMAT (/, 



1 



9H ILLFGAL ,A5,29H. GRAPH WILL BE PLOTTED WITH ,A5, 
5H = 8. ,/ ) 



11 



12 

13 

14 

15 



JWIDE = 8 
GO TO 14 

I F ( IWIDE - 9)13,13,10 
JWIDE = IWIDE 
I F ( IHIGH) 15,16,17 
ITIT = 5HIHIGH 
PRINT 102, ITIT, I T I T 

16 JHIGH = P 
GO TO 19 

17 I F( IHIGH - 15)18,18,15 

18 JHIGH = IHIGH 

19 NODFXAX = MODEXAX 

I F ( MODE X A X )20,27,21 



000918^0 
nonn 1810 
^0001 820 
^ 009 1 8 30 
00001840 
00001850 
00001860 
00001870 
0C001880 
00001890 
00001900 
00001910 
0C001920 
00001930 
00001940 
00001950 
C0001960 
00001970 
00001980 
00001990 
00002000 
ABSOLUTE MAGN00002 010 
00C02020 
0C002030 
n 0002 040 
00002050 
00002060 
) 00002070 

00002080 
00002090 
00002100 
00002110 
00002120 
00002130 
00002140 
00002150 
00002160 
00002170 
00002180 
00002190 
00002200 
00002210 
00002220 
00002230 
00002240 
00002250 
00C02260 
00002270 
00O02280 
0000229 0 
00002300 
0 0002 3 10 
00002320 
00002330 
00002340 



58 



2 r ' 


IT I T= 8HMQDPX&X. 






~99C2359 




PS I NT 1C4* I T I T * IX 






00002360 


104 


FORMAT (/* 9M ILLEGAL *A6* 52 w 


GRAPH a ill 3* PLOTTED WITH 


MODE * 


"0002 3 7 C 




A 1 * 7H£X = r , */) 






"'POO? "3 3 O 




NODEXAX = n 






~ r\ 2 Zi r* 




r , o TO 2 7 






9 C r " 2 4 0 n 


?! 


IF(V,0OFXAX - 2)27*.? 2*22 






P 000 241'' 


22 


I c ( IX'.'P - Jh I Gw )24*24*?? 






00002420 


23 


ITIT = 8HIXUP. 






0 OC 02 4 30 




D ° I N T 1 r 4 * ITIT* IX 






00002440 




NODEXAX = 0 






00002450 




GC TC 27 






00002460 


24 


IF ( IXU P )23*26*26 






*90002470 


26 


JXUP = IXUP 






00002480 


27 


NODE Y AX = MODE Y AX 






O0002490 




I f ( MODE YA X)28*35*29 






00002509 


2 p 


IT IT=8HMODEYAX. 






^ 000 2 5 10 




Dr? I NT 194* ITIT* IY 






0CCP2520 




noqfyax = 0 






0 COO 2 5 30 




GO TO 35 






00002540 


29 


I F ( MODE Y AX - 2)35*30*28 






00002550 


39 


IF(IYPIGHT - JWIOF)32*32*31 






00002560 


3' 


ITIT = 8HIYRIGHT. 






00002570 




PRINT 104* ITIT* IY 






00002580 




NODEYAX = C 






00002590 




GO TO 35 






00002600 


32 


I F ( IYRIGHT) 31*34*34 






00002610 


34 


JYPIGHT = IYRIGHT 






00002620 




INITIALIZE PRIOR TO SCALING 


AND AXIS LOCATING. 




00002630 




IFLAG = 0 FOR PASS WITH XDATA. I FL AG = 1 FOR PASS WITH 


YDATA. 


00002640 


? 5 


DO 2235 IOTA=l*l? 






00002650 


235 


JJTITLE ( IOTA ) = ITITLEf IOTA) 






00002660 




IFLAG = 0 






0CC02670 




SFTA = 0. 






00002680 




SCALE = EXSCALE 






0 0002690 




IAXIS = JYRIGHT 






00002700 




MODE = NODEYAX 






00002710 




I SIZE = JWIDE 






00002720 




I X Y = IX 






00002730 




IYX = IY 






00002740 




AM4X =• AM AX X 






00002750 




AMIN = AMINX 






00002760 




GO TO 52 






00002770 


5n 


IFLAG = 1 






00002780 




BETA = 0. 






00002790 




SCALE = YSCALE 






00002800 




IAXIS = JXUP 






00002810 




MODE = NODEXAX 






00002820 




ISIZE = JHIGH 






00002830 




AMAX = AMAXY 






00002840 




AMIN = AM I NY 






00002850 




I X Y = IY 






00002860 




I vx = IX 






00002870 




CHECK SCALE AND GO TO FIXFD 


OR AUTO SCALE ROUTINES. 




00002880 


5? 


IF(SCALF)53*59*56 






00002890 


53 


PRINT 114* I X Y * IX Y 






00002900 



59 



1 1 4 


F C 9 ~ T (/♦ 9- ILLEGAL * ,1 , ? a ; l*. - 


L - 1 7 •_ ; A I T - 4 




1 


^0 *Al f 7 H-r^ALr, */) 








OO TQ 5Q 




o " .0 0 p G "3 > 




r/Dpccc c | y c r c r a j_ ct pi - ,» t ▼ u ON- - r r.' - c 1 r. •/ 


\~: r *ncf • 




R 4 


rc fS rAL' - 






R 7 


I c ( S r * L ~ 




2^)4^ 


R C 


FALL r C-'-L r !T(S"AL"tI c ra_! 




■> 




^ r A L * - c a ^ T D 0 ■*■ ' ^ • * *■ I 9 C a : 1 " 




", 0 'v O P r Q ~ 




CHF r < AND '*0 V R , JT C AXIS LC r ATjON ! c N* 1 r ~ S c . aR v # cr 


V r ^ c r $ *_ _ 


r> r> -\ j 




T TAG = ^ I- ON ! 9 I A G C A~ L| '.» 1 J c JT 1C, 


JCC c :p . 


r\ r> * r\ 




r^iMCRF - 1 ) ] ^3 ? * 1 n 3 1 * 1 0 - 




^ a r 2 ''i 0 


10 30 


ITAG =0 




? 0 0 0 3 ? 2 0 




GO TO 20^ 




7 o^o 3 0 30 


1031 


PRINT 1 1 r 4 , JYX* IXY, ivy 




'A O ^ -5 O r' 


1104 


FORMAT (/, 5H vo^E* Al »24HAX ^JST NOT 3E 1 UNLESS ,A1 


* 5 7HSCALE IS 


'■.■''no'VJ ^ r 


1 


( A!JTO-S r ALF) • OR*pu MLL 8 C PLOTTED WITH A l j T C * A l , 


7 M - S G 4 L c • , / ) 






GO TQ RQ 




^ ^ •> -j r* -» « 


1^32 


I r ( a P F C ( X V A X - AM T N > - 1 • 0 C -97 ) 1 0 8 3 , 1 r 3 3 , X 03 8 




f\ 0 ^ ^ 3 r> fl n 


1033 


I c- ( A RFC ( AM AY 5 - ! # r ' £: -C7)^H,l^^o,in?.c 




0 00^3 09 n 


1039 


1 r ( A VAX ) 1036* 1" 04, r^T 




0 C^C 3 100 


1034 


I AX IS = I o t Z F / 2 




^0003 1 10 




GO TO 1030 




00003 1 20 


1036 


I AX I S = I SIZE 




00003130 




GO TO 1030 




00003140 


1037 


I AX IS = 0 




00003150 




GO TO 1030 




0000 3 16 3 


1038 


IF (SIGN F( l. , A MAX! - SIGNF( 1 • * A M I N ) 1 1040,1030*1340 




00003170 


1040 


AS I Z F = I S I Z F 




CQ003 1 30 




I A X I S = - AM I N / ( AM A X - A M I N ) * A S I Z F +0,5 




003190 




GO TO 103 0 




0 0 0 ^ 3 2 0 0 




AIJTO SCALE ROUTINE. 




^00^-3710 


59 


IF (MODE - 1) 60 *64, 69 




OOQ03 220 


60 


AMAX = M A X 1 F ( 0 # 9 AMAX) 




00003230 




AMIN = MIN1F (0* * AMIN) 




00003240 


64 


IF(ABSF(AMAX - AMIN) - 1 . OE-97 } 65 , 68 , 68 




00003250 


65 


PRINT 116* IXY, I X Y , IYX 




00C03260 


1 16 


FORMAT (/* 5H ALL *A1*47H VALUES EQUAL. AUTO SCALE 


P0S$H3L C ONLY 


00003270 


1 I F THE *A1,?9H VALUES ARE NON-ZERO AND V0DE.A1* 7HAX 


= 2 . ) 


00003280 




GO TO 1005 




00003290 


68 


ASIZF = I SIZE 




00003300 




SCALE = (AMAX - AM I N ) / AS I Z F 




00003310 




GO TO 83 




00° 03 320 


69 


IF(ARSF(AMAX - AMIN) - 1. OE-97) 70*74*74 




00003330 


70 


IF ( ARSF ( AMAX ) - 1 . OF -9 7 ) 7 1 , 74 , 7 4 




00003 340 


71 


PRINT 118* IXY 




00003350 


118 


FORMAT </, 5H ALL . A 1 . 3 HH VALUES ZERO. AUTO SCALE NOT POSSIBLE. 


) 00003 360 




GO TO 1005 




00003370 


74 


I F ( AMA X - 1 . OE-97) 76* 75* 75 




00003380 


75 


IF( ISIZF - I A X I S ) ? 7 , 76* 71 




00003390 


76 


SCALF1 = 0. 




00003400 




GO TO 7 P 




D0003410 


77 


AXIS * I A X I S 




G 0 0 0 3 4 2 0 




aftze = I^I? C 




00003430 




SCALF1 = AVAX/(ASI/F - AXIS) 




00rn3440 


78 


IF (AMIN + | # or- -97) 79*79*80 




00003450 


70 


I F ( TAXIS)R1, 80*81 




00003450 



60 



nn n n 



80 SCAL c 2 = 

GO TO 62 

81 AXIS = I A x I S 



00003490 

C0003490 



SCALE2 = -AMIN/AXIS 

82 I F ( SC M_ p 1 + SCAL c 2)1994, 1932* 1084 

1982 PRINT 1983* IYX* IYX 

1983 FORMAT (/* 56H NONE 0 f THE PLOT LIES ON 

1IFR *A1*47H — AXIS LOCATION* GRAPH W ILL 3 C 
? 7HAX = *. ,/) 

MO OF = 0 
GO TO 6C 

1984 SCALF = MAX1F ( SCALF 1 * SCALF2) 

83 CALL SCAL^I T ( SCALF , I SC AL1 0 , F a r TOR * 3 ) 



000035C0 
0000 35 1 C 
n 0CC35 20 

THE GRAPH WITH THIS SPEC ! F0000353C 
PLOTTED WITH v 0DE * A 1 * *0093540 

000O35F0 

00003560 

00003570 

00003580 

00993 c 90 



I F ( FACTOR - 5*05)95*65*84 

84 FACTOR = 1 

ISCAL10 = ISCAL10 + 1 
GO TO 90 

85 IF (FACTOR - 2 • 02 ) 87 * 9 7 * 86 

86 FACTOR = 5 
GO TO 90 

87 I F ( FACTOR - 1.01)89,89*89 

88 FACTOR = 2 
GO TO 90 

89 FACTOR = 1 

90 SCALF = FACTOR*10.**ISCAL10 

C COMPUTE AXIS LOCATION IF NECESSARY. AUTO SCALE CASE. 

I F ( MODE - 1)92*91*93 

91 I F ( S 1 GN F ( 1 • * AM A X ) - S I GNF ( 1 . * AM I N ) ) 92 * 94 * 92 

92 IAXIS = -AMIN/SCALE + 0.5 

93 ITAG = 0 
GO TO 203 

94 IF (AMAX )95*95*200 

95 IAXIS = I S I 2 E 
RFTA = -AMAX/SCALF 

I F ( BE T A - l.E+12 >99,99*96 



00003600 

0C003610 

00003620 

00003630 

00003640 

00003650 

00*^3 

C000367Q 

00003680 

00003690 

00003700 

00003710 

00003720 

00003730 

00003740 

00003750 

00003760 

00003770 

00003780 

00003790 

00003800 

00003810 



96 PRINT 120, IXY 

12C FORMAT (/, 15H THF ORIGIN OF ,A1, 43H CANNOT BE OFFSET MORE THAN 
l.OF+12 INCHES. ) 

GO TO 1005 

99 IRFTA = BETA 4 0.5 
PFT A = -IRFTA 

BETA IS THE N'JMBFR OF INCHES OF ORIGIN 
TRU C ORIGIN IS R c L OW OP TO L p FT OF TH C 



SUPPRESS! ON * POS 1 T I VE 
GRAPH. 



00003820 
100003830 
00003840 
00003850 
00003860 
00003370 
I F 00003880 
000C3890 



IF ( BETA + 1. ) 97,97*93 

97 ITAG = 1 
GO TO 203 

200 IAXIS = 0 

BFT A = AM I N / SC AL C 
I F ( BFT A - l.E + 12) 201 *201,96 

201 IBETA = BFT A + 0.5 
BETA = I BFT A 

I F ( BET A - 1. >93 *202,202 

202 ITAG = 1 

RELEASE RESULTS TO REMAINING PART OF PROGRAM. START SECOND 
PASS FOR Y VALUES IF NOT Y F T COMPUTED. 

2 03 I F ( I FLAG ) 2 n 5*2D4,205 



00003900 

00003910 

00003920 

00003930 

00003940 

00003950 

00003960 

00003970 

00003980 

00003990 

00004000 

00004010 

00004020 



61 



c\ n 



204 



205 



r 



206 



9206 

207 

208 

2 0 Q 
210 



211 

212 

213 



214 

MS 



216 

C 

C 



^ALEX = SCAL? 

IX- An = Fan OR 
I X S C 1 " = I 5C A L 1 0 
IXAXIS = I A X I S 
ITAGX = ITAG 
IMZEX = ISIZE 
8ETAX = BETA 
GO TO 50 

T A Y = BETA 
SCALFY = SCALE 
IYFACT = FACTOR 
I v S r 1 o = I C G A L 1 ° 

I Y A X I S = I A X I S 

NOW W 9 I T p r?Ffnc?r)C # 
ITAGY = ITAG 



ISIZEY = 
THIS 
TAPE 
JXT IT ( 1 ) 
JXT I T ( 2 ) 
JXT IT( 3 ) 
JXT IT ( 4 ) 
JX T I T ( 5 ) 
JYTIT ( 1 ) 
JYTIT! 2 ) 
JYT I T ( 3 ) 
JYT IT (4) 
JYTIT! 5 ) 
DC 9206 
JXT I T ( I ) 
J Y T I T ( I ) 



I S I Z F 

COMPLETES CALCULATUN CF SCALE FACTORS ETC. 
RECORDS. FIRST* THE SCALE FACTOR TITLES. 

= 8 H 1 X- 



NOW GENERATE 



= 8H5CALF = 

= I COD c ( SCAL p X ) 

= QH UNITS/ I 
= 8HNCH. 

= 8H1 Y- 

= RHSCALF = 

= ICODF(SCALFY) 

= 8H UN ITS/ I 
= 8HNCH. 

1 = 6,11 
= 8H 
= 8H 

I F ( ITAGX )211*211*207 
IF (BETAX) 208*208*209 
JXT IT ( ) = 8H ADD - 
GO TO 210 

JX T I T ( 7 ) = AH ADD + 

JX T I T ( 8 ) = I CODS ( 9ETAX*SCALFX ) 

JXT I T ( 9 ) = 8H UNITS T 
JX T I T ( 10) = 8HO ALL X 
JXT I T ( 11)= 8HVALUES. 

I F ( ITAGY)216*216*212 
IF(BETAY)213,213,214 
JYT IT ( 7 ) = 8H ADD - 
GO TO 215 



8H ADD + 

ICODE ( 3ETAY*SC4LEY ) 
8 H UNITS T 



JYT I T ( 7 ) = 

JY T I T ( 8 ) = 

JY T I T ( 9 ) = 

JY T I T ( 1 0 ) = 8HO ALL Y 
JYTIT ( 11 )= 8HVALUFS. 

I FONT ( 1 ) = ICONTRL + 4 

INSERT TITLF SIZ^ (02B) AHPAD 
CALL ISHIFT6 (ITITLE* LTI1LF) 

TEST FOR ALL BLANK TITLES. 
ICHECK = 8H 
DO 9075 1=1*6 



OF MAIN TITLE RECORD « 



I F ( ITITLE ( I )- ICH p CK) 9074,9075*9074 



00034030 

CC004 

00"C4 050 
0000 a 060 
OCOO^MO 
3CO04 080 
00004090 
C0004100 
00004110 
MC04 120 
00004130 
00004 140 
00004150 
00004160 
0 0004 1 70 
0000418" 
00004190 
00004200 
C 0004 2 10 
00004220 
00004230 
00004240 
00004250 
00004260 
00004270 
00004280 
00004290 
00004300 
00004310 
00004320 
00004330 
00004340 
00004350 
00004360 
00004370 
00004380 
00004390 
00004400 
00004410 
00004420 
00004430 
00004440 
00004450 
00004460 
00004470 
00004480 
00004490 
00004500 
00004510 
00004520 
00004530 
00004540 
00004550 
00004560 
00004570 
00004580 



62 



9074 


IF ( ITT TLP C I » ) 9080,9075,9^80 


90004590 


9075 


CONTINUE 


C00046CC 




IT1 = 1 


90004610 




I CONT ( 1 ) = I OONT ( 1 ) - 1 


00904620 




GO TO 9081 


00004630 


9080 


IT 1 = 0 


00004640 


9081 


DO 9085 1=7*12 


00004650 




IF (ITITLE(I) - ICHECK) 9084,9085,9084 


00004660 


9084 


IF ( ITITLE( I ) )9090, 9085, 9090 


00004670 


9085 


CONTINUE 


00004680 




I T 2 = 1 


00004690 




I CONT ( 1 ) = I CONT ( 1 ) - 1 


00004700 




GO TO 9091 


00004710 


9090 


IT2 = 0 


00004720 




NOW gfneratf axes records. 


00004730 


909] 


LFTMGN = 0*100 


00004740 




IBOTMGN = 1*100 


00004750 




I 1 - = LFTf.GN 


00004760 




jH = IBOTMGN + I YAX I S* 1 00 


00004770 




LH = I S IZEX*100 


00004780 




IHL = LFTMGN + ISIZEX*100 - 107 


00004790 




KAXIS(l) = IPACK12(IH,JH,LH,IHL) 


00004800 




JHL = JH - 13 


00004810 




I HL2 = -100 


00004820 




IVH = (ISIZEX -IXAXIS - 1 ) * I XFACT 


00004830 




I VH2 = -IXFACT 


00004840 




KAX I S ( 2 ) = IPACK12(JHL,IHL2»IVH,IVH2) 


00004850 




NH = ISIZFX 


00004860 




ISH = 8H 14 


00004870 




IV = LFTMGN + IXAXIS*100 


00004880 




JV = IBOTMGN 


00004890 




KAX I S ( 3 ) = IPACK12INH, ISH, IV, JV) 


00004900 




LV = I S IZEY*100 


00004910 




IVL = IV - 3 


00004920 




JVL = IBOTMGN + ISIZEY*100 - 107 


00004930 




JVL2 = -100 


00004940 




KAX I S ( 4 ) = IPACK12(LV, I VL » JVL , JVL2 ) 


00004950 




IVV = (ISIZEY - I YAX IS - 1)*IYFACT 


00004960 




I VV2 = -IYFACT 


00004970 




INV = ISIZEY 


00004980 




ISV = 8H 11 


00004990 




K AX I S ( 5 ) = IPACK12(IVV,IVV2,INV*ISV) 


00005000 




NOW GENERATE CURVES# 


00005010 




SCX = 100 • /SCALEX 


00005020 




SCY = 100# /SCALEY 


00005030 




EX AX I S= I XAX I S* 1 0 0 


00005040 




YAXIS = I YA X I S* 100 


00005050 




ALFTMGN = LFTMGN 


00005060 




BOTMGN = IBOTMGN 


00005070 




SHIFTX sEXAXIS - BFTAX*100# + ALFTMGN 


00005080 




SHIFTY = YAXIS - BFT AY*100 • + ROTMGN 


00005090 




EXSIZF* ISIZEX* 100 + LFTMGN + 60 


00005100 




SIZEX * LFTMGN - 60 


00005110 




YSIZE = I S IZEY*100 + IBOTMGN + 70 


00005120 




SIZEY = IBOTMGN - 70 


00005130 




ICURV(l) = 0 


00005140 



63 



240 


I F { I POINT ) 90 10, 900 7, 9010 


00005150 


9007 


IF (XMODF(NUMPTS»2 ) ) 9 70 0 * 9 70 1 * 9 700 


00005160 


9700 


ISWITCH = 1 


00005170 




60 TO 242 


00005180 


9701 


ISWITCH = 2 


00005190 


242 


INUM = (NUMPTS + 1 ) /2 


00005200 




DO 244 1=1*1 NUM 


00005210 




Cl = X ( 2* I - 1 ) *scx + shiftx 


00005220 




C2 = Y < 2*1 — 1) *SCY + SHIFTY 


00005230 




IF ( I -I NUM ) 241 *9241 ,241 


00005240 


9241 


GO TO (9242,241) , ISWITCH 


00005250 


9242 


C3 = Cl 


00005260 




C4 = C2 


00005270 




60 TO 9243 


00005280 


241 


C3 = X ( 2* I ) *SCX + SHIFTX 


00005290 




C4 = Y( 2*1 ) *SCY + SHIFTY 


00005300 


9243 


Cl = MIN1F(C1*EXSIZE) 


00005310 




IC1= MAX 1 F ( C 1 , SIZEX) 


00005320 




C2 = MIN1F(C2* YSIZE) 


00005330 




IC2 = MAXI F ( C2 # SIZEY) 


00005340 




C3 = MIN1F ( C3*EXSIZE ) 


00005350 




I C3 = MAXI F { C3 * SIZEX) 


00005360 




C4 = M I N1F ( C4 * YSIZE) 


00005370 




I C4= MAXI F ( C4 * SIZEY) 


00005380 


244 


ICURVU + 1) = IPACK12 ( IC1* IC2* IC3* IC4) 


00005390 




I I = INUM + 3 


00005400 


246 


CALL I PACKL 1 ( LABEL * LABEL 1 * IDUMMY) 


00005410 




I CURV ( II-l) = LABEL1 


00005420 




ICURV(II) = ICURV4 


00005430 


9010 


I F ( MODCURV - 1)247*247*9015 


00005440 


247 


CALL IREADY (IDUMMY) 


00005450 




I F ( IDUMMY) 5000, 1260,5000 


00005460 


1260 


CALL IWRITE (ICONT* IDUMMY, 1) 


00005470 




I F ( IDUMMY) 5000,260*5000 


00005480 


260 


CALL IWRITE (JXTIT, IDUMMY, 11) 


00005490 




I F ( IDUMMY) 5000,261 ,5000 


00005500 


261 


CALL IWRITE (JYTIT, IDUMMY, 11) 


00005510 




I F ( IDUMMY) 5 000, 265, 5000 


00005520 


265 


I F ( IT1 )9269, 9268, 9269 


00005530 


9268 


CALL IWRITE (LTITLF, IDUMMY, 7) 


00005540 




I F ( IDUMMY)50C0, 9269, 5000 


00005550 


9 2 b 9 


I F ( IT2 )9271, 9270, 9271 


00005560 


9270 


CALL IWRITE ( LT I TLE ( 7 ) , IDUMMY, 7) 


00005570 




IF ( IDUMMY) 5 000, 92 7 1,5 000 


00005580 


9271 


CALL IWRITE (KAXIS, IDUMMY, 5) 


00005590 




I F ( IDUMMYJ5000, 9015, 5000 


00005600 


9015 


I F ( I PO I NT ) 9 020, 270, 902 0 


00005610 


270 


CALL IWRITE (ICURV, IDUMMY, II) 


00005620 




I F ( IDUMMY) 5000,9020,5000 


00005630 


9020 


I F ( MODCURV - 1)272,272,9025 


00005640 


272 


IFUGRID - 1)9025,273,9025 


00005650 




GENERATE GRID IF CALLED FOR* 


00005660 


273 


1X100 = I S I ZEX*100 


00005670 




I Y 100 = ISIZEY*ino 


00005680 




NEXT 1 = IBOTMGN 


00005690 




NFXT2 = LFTM6N + 1X100 


00005700 



64 



JGR I D ( 1 ) = 0 
DO 1274 J=1 *19*2 
JGR I D ( J+ 1 ) = I PACK 12 
I F ( NEXT 1 - IBOTMGN - 

1273 NEXT1 = NEXT 1 + 100 
JGR I D ( J+2 ) = I PACK 12 
I F ( NEXT 1 - IBOTMGN - 

1274 NEXT 1 = NEXT 1 + ICO 

1275 JGR I D ( J+2 ) = I PACK 12 

1276 JGR I D ( J+3 ) = ICURV3 
JGR I D ( J+4 ) = ICURV4 

CALL IWRITE( JGRIDt I DUMMY ♦ 
I F ( I DUMMY) 5000 *1277 *5000 

1277 NEXT 1 = LFTMGN 

NEXT2 = IBOTMGN + IY100 
DO 1279 J=l*ll*2 



(LFTMGN* NEXT 1 » NEXT2 # NEXT 1 ) 
IY100 ) 1273*1275* 1275 

( NEXT 2 ♦ NEXT 1 • LFTMGN* NEXT 1 ) 
I Y100 ) 1274*1276 *1276 

(NEXT2* NEXT1* NEXT2* NEXT 1 ) 



J+4) 



I F ( I DUMMY ) 5 000 *902 5* 50 00 
9025 IF ( IPOINT )9030*276*9030 
I GENERATE POINT PLOT RECORDS 

I OUT = 0 

CALL IPACKL1 (LABEL* LABEL 1 ♦ LABEL2 ) 
DO 9050 1=1 *NUMPTS 
Cl = X( I)*SCX + SHIFTX 
C2 = Y( I ) *SCY + SHIFTY 
I F ( C 1 - EXSIZE)9031*9031*9034 
I F ( C2 - YSIZE)9032*9032*9034 
I F ( C 1 - SIZFX)9034*9033*9033 
I F ( C2 - SIZEY )9034*9035*9035 
I OUT = I OUT +1 
GO TO 9050 
IC1 = Cl 



IF CALLED FOR. 



9030 



9031 

9032 

9033 

9034 

9035 



00005710 

000C5720 

00005730 

00005740 

00005750 

00005760 

00005770 

00005780 

00005790 

00005800 

00005810 

00005820 

00005830 

00005840 

00005850 

00005860 





JGRIDt J+l) 


= I PACK 12 


( NEXT 1 1 


IBOTMGN* NEXT1 


* NEXT 2 ) 


00005870 




I F ( NEXT 1 - 


LFTMGN - 


1X100) 1278»1280*1280 




00005880 


1278 


NEXT 1 = NEXT 1 + 100 








00005890 




JGRIDt J+2) 


= I PACK 12 


t NEXT 1 * 


NFXT2* NEXT1. 


IBOTMGN) 


00005900 




I F ( NEXT 1 - 


LFTMGN - 


1X100)1279*1281*1281 




00005910 


1279 


NEXT 1 = NEXT 1 + 100 








00005920 


1280 


JGRIDt J+2) 


= I PACK 12 


( NEXT 1 » 


NFXT2, NEXT 1 « 


NEXT2) 


00005930 


1281 


JGRIDt J+3) 


= ICURV3 








00005940 




JGRIDt J+4) 


= ICURV4 








00005950 




CALL IWRITE 


( JGRID. 


IDUMMY* 


J+4) 




00005960 



00005970 

00005980 

00005990 

00006000 

00006010 

00006020 

00006030 

00006040 

00006050 

00006060 

00006070 

00006080 

00006090 

00006100 

00006110 





GO TO (9036 


*9037*9038*9039 


*9040) • 


IPOINT 




00006130 


C 


GENERATE CROSS. 










00006140 


9036 


ICURVI2) = 


IPACK12 


( I Cl-5 ♦ 


IC2-5* 


IC1 + 5 * 


IC2+5) 


00006150 




I CURV ( 3 ) = 


IPACK12 


( IC1 ♦ 


IC? * 


I Cl-5 ♦ 


+C2+5) 


00006160 




I CURV ( 4 ) = 


IPACK12 


( I Cl +5 * 


IC2-5. 


IC1+5* 


IC2-5 ) 


00006170 




GO TO 9041 












00006180 


C 


GENERATE PLUS. 










00006190 


9037 


I CURV ( 2 ) = 


IPACK12 


( IC1 * 


IC2-5* 


IC1 * 


IC2+5) 


00006200 




I CURV ( 3 ) = 


IPACK12 


t IC1 * 


IC2 * 


I C 1-5 ♦ 


IC2 ) 


00006210 




I CURV ( 4 ) = 


IPACK12 


( IC1+5* 


IC2 * 


IC1+5* 


IC2 ) 


00006220 




GO TO 9041 












00006230 


C 


GENERATE SQUARE 


. 








00006240 


9038 


I CURV ( 2 ) = 


IPACK12 


( IC1+4* 


IC2-4* 


I C 1 +4 * 


IC2+4) 


00006250 




I CURV ( 3 ) = 


IPACK12 


( IC1-4. 


IC2+4* 


I C 1-4 » 


IC2-4) 


00006260 




I CURV ( 4 ) = 


IPACK12 


( IC1+4* 


IC2-4* 


I C 1 + 4 * 


IC2-4) 


00006270 
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GO TO 9041 

C GENERATE DIAMOND. 



9039 


I CURV ( 2 ) = 


IPACK12 


( I C 1 + 5 » 


I C 2 , 


IC1 * 


IC2+5) 




I CURV ( 3 ) = 


IPACK12 


( IC1-5. 


IC2 , 


IC1 , 


I C 2 — 5 ) 




I CURV ( 4 ) = 

GO TO 9041 


IPACK12 


< IC1+5* 


IC2 * 


I C 1 + 5 , 


I C 2 ) 




GENERATE TRIANG 


LE. 








9040 


I CURV ( 2 ) = 


IPACK12 


( IC1+5. 


IC2-3* 


IC1 . 


IC2+6) 




I CURV ( 3 ) = 


I PACK 1 2 


( I C 1-5 * 


I C 2-3 * 


I C 1 + 5 * 


IC2-3) 




I CURV ( 4 ) = 


I CURV ( 3 ) 










9041 


I F ( I - NUMPT S ) 9043 « 


9042 *9043 






904? 


I CURV ( 5 ) = 


LABEL 2 











GO TO 9046 

9043 I F < I - 1)9045*9044*9045 

9044 I CURV ( 5 ) = LABEL 1 



GO TO 9046 

9045 I CURV ( 5 ) = ICURV3 

9046 I CUR V ( 6 ) = I CURV4 

CALL I WRITE (I CURV* I DUMMY , 6) 

I F < I DUMMY )5000*9050*5000 
9050 CONTINUE 

I F ( I OUT ) 9048 *276*9048 
9048 PRINT 9104, IOUT 

9104 FORMAT (/* IX* 12* 29H POINT(S) WERE OFF THE GRAPH. ,/) 
C SET UP RETURN. 

276 IF(MODCURV)277, 278*277 

277 I F ( MODCURV - 3)279*278*279 

278 ITEST = 0 

PRINT 130*( JJTITLEU ) *1 = 1.12) 

130 FORMAT (/. 19H GRAPH TITLED . . . 6A 8 . / . 1 9X * 6A 8 , 

1 24H . . HAS BEEN PLOTTED. ./.1H0) 

I DUMMY = ITYP2 ( I DUMMY ) 

I F ( I DUMMY ) 5 670.656.5670 
656 LAST = 0 
RETURN 

279 I TEST = 1 

IDUMMY = I CLOCK ( I DUMMY ) 

LAST = 0 
RETURN 

THESE ARE THE NORMAL RETURNS. 

NOW SET UP THE RETURN FOLLOWING A TAPE ERROR. 

5000 I F ( MODCURV - 1)5001*5001*5002 

5001 IDUMMY = I TYPE 1 ( I DUMMY ) 

GO TO 247 

5002 PRINT 5100 

5100 FORMAT (/* 36H TAPE ERROR IN WRITING GRAPH OUTPUT. ) 
IDUMMY = ITYPE1 (IDUMMY ) 

GO TO 1007 

5670 IDUMMY = I TYPEl ( I DUMMY ) 

END 



SUBROUTINE IREADY mhiimmY) 

SELECTS TAPE 8 (. LL LOOP UNTIL READY). WRITES EOF ON 8. 
MACHINE LANGUAGE WILL NOT BE NECESSARY IN FORTRAN 62-3. 
LOCUFIVE = 5 ) • 



00006280 

00006290 

0000630C 

00006310 

00006320 

00006330 

00006340 

00006350 

00U06360 

00006370 

00006380 

00006390 

00006400 

00006410 

00006420 

00006430 

00006440 

00006450 

00006460 

00006470 

00006480 

00006490 

00006500 

00006510 

00006520 

00006530 

00006540 

00006550 

00006560 

00006570 

00006580 

00006590 

00006600 

00006610 

00006620 

00006630 

00006640 

00006650 

00006660 

00006670 

00006680 

00006690 

00006700 

00006710 

00006720 

00006730 

00006740 

00006750 

00006760 

00006770 

00006780 

00006790 

00006800 

00006810 

00006820 

00006830 
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FXF (520413) 


EX F 7 (520005). 


SELECT READ AND WAIT TAPE. 


1NEX 


FXF7(0005^P) 


SLJ (1RDY). 


E v i t ON r H 3 ACTIVF. 




EXF7 ( 529C0B ) 


SLJ (1NEX). 


EXIT ON TAPF READY. 




LOA ( I F I VF ) 


SVJ (1RUF). 


t-rminatf 


13U“ 


EXF 5 ( N ) , 




BUFFER. 


1 RDY 


EXF ( 52041 B ) 
ENA (0). 


E XF 7 ( 52 00 OB ) , 


SELECT AND WAIT TAPE 8. 
clear A. 




EXF (020003) 


EXF (520063). 


STOP CLOCK AND BACKSPACE 8. 




EXF7( 52001B) 
-EXF7 ( 520003 ) , 


SLJ ( TEND ) . 


EXIT IF NOT AT LOAD POINT. 
WAIT TAPE 8. 




EXF7( 520073) 


SLJ ( 1 EOF ) . 


EXIT IF NO EOF. 




ENA (IDUMMY) 


SAU ( 2BUF ) . 


MOVE 




INA (1) 


STA (IFIVE). 


FORWARD 


2BUF 


EXF5 ( N ) 


E XF 7 ( 5200CB ) . 


OVER RECORD. 


1 EOF 


ENA (0) 


EXF7( 00C61B). 


clear a. wait ch 6. 




EXF ( 6204 1 B ) 


EXF7 ( 620003) . 


SELECT WRITE AND WAIT TAPE. 




FXF ( 62041 B ) 


EXF7 ( 62000B) . 


SELECT AND WAIT TAPE 8. 




FXF (62003B) 


EXF7 ( 620003 ) . 


WRITE EOF AND WAIT. 


1 END 


EXF7 ( 62007B ) 
STA ( IDUMMY ) • 


ENA( 10) . 


EXIT ON NO END OF TAPE. 



END 



SUBROUTINE I WR I TE ( I S TART * IDUMMY, IWORDS) 

WRITE RECORD OF IWORDS, STARTING WITH ISTART. PUT IDUMMY = 0 
IF RECORD CORRECTLY WRITTEN* OTHERWISE SET NON-ZERO, 

MACHINE LANGUAGE WILL NOT RE NECESSARY IN FORTRAN 62-3. 





LOC ( I S I X = 6) . 
-F.X F7 ( 00061 B ) . 
EXF ( 6 204 IB ) 
FXF ( 6 204 IB ) 
ENQ ( 1110) • 


EXF7 ( 6200QB ) . 
EXF7 ( 62000B ) • 


WAIT CH 6. 

SELECT WRITE* WAIT TAPE. 
SELECT AND WAIT TAPE 8. 
SET COUNTER. 


1AGN 


ENA ( ISTART ) 
SAL ( 1 BUF ) 


I NA t 1 ) . 

ADD ( IWORDS) . 


STARTING ADDRESS. 
TERMINAL ADDRESS. 


1 BUF 


ST A ( I S I X ) 

ENA ( 0 ) 

EXF7 ( 62007B ) 
FXF7 ( 6 20O3B ) 
EXF7( 620043) 


EXF6I N I . 

EXF7 ( 62000B ) . 
SLJ ( 1END) • 
SLJ ( 2 AGN ) • 
SLJ t 2 END ) • 


BUFFER OUT. 

CLEAR A. WAIT TAPE 8. 
EXIT IF NO END OF TAPE. 
EXIT IF NO PARITY ERROR. 
EXIT IF LENGTH ERROR. 


2 AGN 


FXF ( 620063 ) 
ORS (3) 


EXF7 ( 62Q00B ) • 
Q JP 1 ( 1 AGN I . 


BACKSPACE AND WAIT. 
TRY WRITE 3 TIMES. 


1 END 
2END 


FXF (620O3B) 
STA (IDUMMY). 


ENA (10). 


WRITE EOF, NON ZERO A. 
STORE RESPONSE. 



END 



FUNCTION ITYP2 (IDUMMY) 

TYPE WORD GRAPH, 

WILL NEED REWRITING IN FORTRAN 62-3, 

CON ( LC = 5 7B * Ml = 4 5 1 3 1 23 01 50 5 000 OB )• 

LOC ( I TWO = 2 ) • 

-EXF7 ( 00061B ) , WAIT CH 6. 



FXF (62041R) 
FXF (62041ft) 
FXF (62003B) 
ENA (0) 

EXF7( 620076) 
STA (ITYP2). 



F XF 7 ( 62000B) , 
FXF 7 ( 62O00B), 
EXF (010003), 
EXF7 ( 62000B ) . 
ENA (10), 



SELECT AND WAIT TAPF, 
SELECT AND WAIT TAPE 8, 
WRITE EOF, START CLOCK, 
CLEAR A, WAIT TAPE 8. 
EXIT IF NO END OF TAPE, 
STORE RESPONSE. 



0 0006840 
00006850 
O0006860 
00006870 
00006880 
00006890 
00006900 
00006910 
00006920 
00006930 
00006940 
00006950 
00006960 
00006970 
00006980 
00006990 
00007000 
00007010 
00007020 
00007030 
00007040 
00007050 
00007060 
00007070 
00007080 
00007090 
00007100 
00007110 
00007120 
00007130 
00007140 
00007150 
00007160 
00007170 
00007180 
00007190 
00007200 
00007210 
00007220 
00007230 
00007240 
00007250 
00007260 
00007270 
00007280 
00007290 
00007300 
00007310 
00007320 
00007330 
00007340 
00007350 
00007360 
00007370 
00007380 
00007390 
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-EXF 7 ( 0002 IB ) . 




WAIT CH 2. 


00007400 


EXF7I 11141B) 


SLJ (1TYP). 


EXIT IF UPPER CASE. 


00007410 


EXF (211O0B) 


ENA ( LC+1 ) • 


TYPE 


00007420 


STA ( I TWO) 


EXF2( LC) . 


LOWER CASE. 


00007430 


-EXF7I0C021B). 




WAIT CH 2. 


00007440 


1 T Y D FXF (211008) 


ENA (Ml+1). 


TYPE 


00007450 


STA ( I TWO ) 


E X F 2 ( M 1 ) . 


GRAPH 


00007460 


END 






00007470 








00007480 


FUNCTION ITYPE1 


( I DUMMY ) 




00007490 


REWIND TAPE 


8. REQUEST NEW TAPE. AND WAIT TILL READY. 


00007500 


WILL NEED REWRITING IN FORTRAN 62-3* 


00007510 


CON ( LC = 57B. Ml 


= 4515112030242004B. M2 = 1 1033022040620 3 IB . 


00007520 


1 M3 = 0401301520043342B) . 




00007530 


RSV ( ME SS = 3)* 






00007540 


LOC { I TWO =2). 






00007550 


-EXF7(00061B) . 




WAIT CH 6. 


00007560 


EXF (62041B) 


EXF 7 ( 62000B ) . 


SELECT AND WAIT TAPE. 


00007570 


EXF { 6 204 1 B ) 


EXF7 ( 62000B ) . 


SELECT AND WAIT TAPE 8. 


00007580 


EXF ( 62003 B ) 


EXF7 ( 620006 ) « 


WRITE EOF AND WAIT. 


00007590 


EXF ( 62007B ) 


EXF7 ( 0002 IB ) # 


REWIND WITH INTERLOCK. WAIT CH 


200007600 


EXF7 ( 11141B) 


SLJ (1TYP). 


EXIT IF UPPER CASE. 


00007610 


EXF (21100B) 


ENA (LC+1). 


TYPE 


00007620 


STA ( I TWO ) 


EXF 2 ( LC ) . 


LOWER CASE. 


00007630 


1TYP LDA (Ml) 


STA (MESS). 




00007640 


LDA (M2) 


STA ( MESS+1 ) . 




00007650 


LDA (M3) 


STA ( MESS+2 ) . 




00007660 


-EXF7 ( 0002 IB ) • 




WAIT CH 2. 


00007670 


EXF (21100B) 


ENA ( MESS+3 ) . 


TYPE 


00007680 


STA ( I TWO ) 


EXF2 ( MESS ) . 


MESSAGE. 


00007690 


-EXF7 ( 0006 IB ) • 




WAIT CH 6. 


00007700 


-EXF7 ( 62000B ) • 




WAIT TAPE. 


00007710 


EXF ( 6204 IB ) 


EXF7 ( 62000B ) . 


SELECT AND WAIT TAPE 8. 


00007720 


EXF ( 0 10 00B ) • 




START CLOCK. 


00007730 


END 






00007740 








00007750 


FUNCTION ICODE (ANUMBER) 




00007760 


CODES ABSOLUTE VALUE OF A FLOATING POINT NUMBER (BETWEEN 


00007770 


1.0E-100 AND 


1.0E+100) INTO 8 


-CHARACTER BCD WORD OF THE FORM 


00007780 


1.23E+45. ICODE = 8H0.00E+00 


IF MAGNITUDE OUT OF RANGE. 


00007790 



CHECK AVAILABILITY OF LIBRARY FUNCTIONS IN FORTRAN 62-3. 00007800 

DIMENSION 11(8) 00007810 

BNUM3ER = ABSF ( ANUMBFR ) 00007820 

I F ( BNUMBER - 1 . 0E + 100 ) 7 .6 .6 00007830 

7 I F (BNUMBER - 1 . 0E-100 ) 6 .6 . 2 00007840 

6 ICODE = 8H0.00E+00 00007850 

RETURN 00007860 

THIS IS ERROR EXIT. 00007870 

2 CALL SCALE I T (BNUMBER. ISCAL10. FACTOR. 3) 00007880 

ISIGNSC = XSIGNFt 1.ISCAL10) 00007890 

ISCAL10 = XABSFf ISCAL10) 00007900 

I FACT = FACTOR* 100.001 00007910 

11(8) = XMODF ( ISCAL10.10) 00007920 

I I ( 7 ) = I SCAL10/ 10 00007930 

I F ( ISIGNSC) 4.3.3 00007940 

3 I I (6) = 8H + 00007950 
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4 

5 



C 



1 



2 

3 



4 

5 



6 

7 

8 
9 

10 

11 



12 

13 

i 



1NFX 



G ) TO 5 

I I (6) = 8H 

11(5) = 8H E 

11(4) = XMODF ( I FACT * 10 ) 

11(3) = (XMODM IFACTflOC) )/10 

11(2) = 8 H 

11(1) = I FACT / 1 0 0 

CALL IPACK (II, IPACKED) 



I CODE = 

RFTUPN 

END 



I PACKED 



SUBROUTINE SCALE I T ( ANUM3EP * 1 SCAL 1 0 ♦ F ACTOR t MODE ) 

FINDS FACTOR (BETWEEN 1.0 AMD 9.99...) AND SCALE OF 10 AS 
DEFINED BY ANUM3ER = F ACTOR* 1C . * * I SC Al 1 0 . 

MODE IS THE NUMBER OF SIGNIFICANT FIGURES REQUIRED. THIS 



:.0007?60 
00007970 
0C0C7980 
00007990 
00008000 
00008010 
0000802C 
0C0C8030 
00008040 
00008 050 
00008060 
00008070 
00008080 
00008090 
00008100 
MUST 0C008110 



BE BETWEEN 1 AND 10 CR IT WILL BE D UT EQUAL TO SIX. 
CHECK AVAILABILITY OF LOGl^F IN FORTRAN 62-3. 
ISCAL10=LOG10F( ANUYBFR) 

FACTOR = ANUMBER/10.**ISCAL1G 
IF(faCTOR - C.l ) 1 * 2 * 2 
FACTOR = FACTOR*10C. 

I SCAL 1 0 = I SCAL 10 - 2 
GO TO 8 

I F ( FACTOR - 1*0)3, 8*4 
FACTOR = F ACT OR* 1 0 • 

I SCAL 1 0 = I SCAL 1 0 - 1 
GO TO 8 

I F ( FACTOR - 100.0)6*5*5 

factor = factor/ioo. 

I SCAL 1 0 = I SCAL 10 + 2 
GO TO 8 

IF(FACT0R - 10.0)8*7*7 
FACTOR = FACT0R/10 . 

ISCAL10 = I SCAL 10 + 1 
IF ( MODE ) 9 * 9 * 1 0 
MODE = 6 
GO TO 11 

I F ( MODE - 1011*11*9 

I FACTOR = FACTOR* 10.** (MODE - 1) + 0.5 

FACTOR = I FACTOR 

FACTOR = FACT0R/10.**( MODE - 1) 

I F ( FACTOR - 10. ) 13 *12*12 
FACTOR = 1* 

I SCAL 1 0 = I SCAL 1 0 + 1 

RETURN 

END 

SUBROUTINE IPACK (II* IPACKED) 

TAKES 8 SIX-BIT WORDS AND PACKS THEM LEFT TO RIGHT 
IN IPACKED. IF WORD IS ZERO* 12B IS SUBSTITUTED. 
CONVERT TO CODAP FOR FORTRAN 62-3. 

CONUZERO = 12B). 

S I U 1 ( I SAVE ) ENI 1(8). 

LDA 1 ( I I ) AJPK2NEX). 

LDA ( I ZERO ) . 



00008120 

00008130 

00008140 

00008150 

00008160 

00008170 

00008180 

00008190 

00008200 

00008210 

00008220 

00008230 

00008240 

00008250 

00008260 

00003270 

00008280 

00008290 

00008300 

00008310 

00008320 

00008330 

00008340 

00008350 

00008360 

00008370 

00008380 

00008390 

00008400 

00008410 

00008420 

00008430 

00008440 

00008450 

00008460 

00008470 

00008480 

00008490 

00008500 

00008510 



69 



* n n * n n * * n 



2NEX 



C 



1NEX 



2NEX 



LRS (6) 


I N I 1 ( -2 ) • 




C 0008 5 20 


ISKl(-l) 


SLJ (1NEX). 




00008530 


STO ( IPACKED) 


LIU1( ISAVE). 




00C08540 


END 






00008530 








00008560 


SUBROUTINE ISHIFT6 (ITITLE. LTITLE) 


00008570 


INSERTS 02B AHEAD OF 6-WORD TITLE RECORD. 


00008580 


WILL HAVE TO BE CONVERTED TO CODAP IN FORTRAN 62-3. 


00008590 


WATCH ARRAY INDEXING IN FORTRAN 62-3. 


00008600 


CON ( I BLANK = 2020202020202 C20R ) . 




00008610 


SIU1( ISAVE ) 


ENI 1 ( 1 ) . 


SAVE INDEX* SET COUNT. 


00008620 


ENA (2). 




ENTER 02B. 


000C8630 


LDOl ( ITITLE) 


LLS (42). 


PERFORM 


00008640 


ST A 1 ( LTITLE) 


LLS (6). 


SHIFTING. 


00008650 


ISK1(6) 


SLJ ( 1 NEX ) • 


CHECK IF COMPLETE. 


00008660 


LDO ( I BLANK ) 


LLS (42). 


COMPLETE LAST WORD. 


00008670 


ENI 1( 7) 


STAKLTITLE). 


STORE LAST. 


00008680 


ENA (2). 




REPEAT 


00008690 


LOOK ITITLE) 


LLS (42). 


FOR 


00008700 


I N 1 1 ( 1 ) 


ST A 1 (LTITLE) • 


SECOND 


00008710 


I N I 1 ( - 1 ) 


LLS (6). 


TITLE 


00008720 


ISK1 ( 12) 


SLJ (2NEX). 


LINE. 


00008730 


LDO ( I BLANK ) 


LLS (42). 




00008740 


ENI1( 14) 


STAKLTITLE). 




00008750 


LI Ul( ISAVE). 




RESTORE INDEX. 


00008760 


END 






00008770 








00008780 


FUNCTION IPACK12 ( I ONE . I 2 . I 3 . I A ) 




00008790 


PACKS FOUR 


12-BIT WORDS INTO 


ONE 48-BIT WORD. 


00008800 


WILL REQUIRE CONVERSION TO CODAP IN FORTRAN 62-3. 


00008810 


LDA ( IONE ) 


LDO ( 12 ) . 




00008820 


QLS( 36) 


LLS( 12 ) . 




00008830 


LDO ( 13) 


OLS ( 36) . 




00008840 


LLS( 12) 


LDQ( 14 ) . 




00008850 


OLS ( 36 ) 


LLS ( 12 ) • 




00008860 


ST A ( I PACK 12). 






00008870 


END 






00008880 








00008890 


SUBROUTINE IPACKL1 (LABEL* LABFL1 


* LABEL2) 


00008900 


packs two 


4-CHARACTER LABFLS. 




00008910 


USE DECODE/ENCODE IN FORTRAN 


62-3. 


00008920 


CONtIFLAG = 37773777B ) . 




00008930 


LLSI24) 


STA ( LABEL 1 ) • 




00008950 


LDA ( I FLAG ) 


LDO (LABEL). 




00008940 


LDA ( I FLAG ) 


LLS( 24) . 




00008960 


STA ( LABEL2 ) . 






00008970 


END 






00008980 








00008990 


FUNCTION I CLOCK ( I DUMMY ) 




00009000 


FXF (01000R). 




START CLOCK. 


00009010 


END 






00009020 


IC2 = C2 






00006120 


END 






00009030 


END 
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APPENDIX B 



PROGRAM NOMENCLATURE 

The following listing of program symbols is provided to help the 
reader understand symbols used in the digital computer program* Symbols 
identified in section 3 as well as symbols readily identifiable are not 
listed in this appendix, 

X(Y) individual strain gage readings 



1(1) 


percentage of depth of individual strain gage readings 


AZ 


the zero depth intercept of strain axis 


ANG(L) 

SMAX(L) 


the angle of Om** from gage leg number 1 

0 rnA* 


SMIN(L) 


OmiN 


THETA(L) 


Op 


SIGXX(LL) 


experimental value of (Tn 


TRUXY(L) 


experimental value of Tky 


PHI(L) 


the angle 0 


XSXX(L) } 
YSXX(L) J 


components of experimental values of <Jn in proper form 
for plotting 


XTXY(L) ' 
XTXY(L) 


) components of experimental values of 'txy in proper form 
1 for plotting 


zx(i,i) 1 
zx(7,l) J 


1 seven coefficients used in computing smoothed value of On 


XY(1,1) 'j 

XY(7,1) J 


| seven coefficients used in computing smoothed value of ^txy 


ZZZ 


the angle in radians on the circumference of the piping between 
two consecutive smoothed values of (Tn or Txy 


PHI 1(1) 
SIGN 7(1) 


the angle in radians from rosette number 12 of a smoothed 
value of GL or'txy 

• 

smoothed value of On . 
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TAUN 7(1) 



smoothed value of'Txy 



XSNR(I) 

YSNR(I) 

XTNR(I) 

YTNR(I) 



components for plotting zero stress circle on On plot 
components for plotting zero stress circle on Tky plot 



XSN7(I) components of smoothed value of Gw in proper form for 
YSN7(I) j plotting 

XTN7(I) *1 components of smoothed value of in proper form for 
YTN7(I) J plotting 
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APPENDIX C 



DESCRIPTION OF SUBROUTINE GAUSS 2 



A. IDENTIFICATION 

■ TITLE: Solution of Simultaneous Linear Algebraic Equations 

CO-OP ID: F2 
FORTRAN: Gauss 2 

CATEGORY: Simultaneous Linear Equations 

PROGRAMMER: C. B. Bailey and Mary Haynes 

DATE OF COMPLETION: April 1963 

B. RJRPOSE: To solve one or more sets of linear algebraic equations using 

Gaussian elimination with row pivoting and back substitution. 

C. USAGE: 

1. Calling sequence 

Call Gauss 2 (N, M, EP, A, B, X, KER) 

N - order of the matrix in the equation Ax = b, that is, the 
number of linear equations (maximum = 50). 

M- - number of vectors b for which solution vectors x are to be 
obtained, that is, the number of sets of linear equations (Maximum = 60). 



EP - matrix condition parameter (see Mathematical Method). 



A - the elements of the matrix of coefficients of the equations. 



stored in the fora 




i = 1, 2 



, ...» 



N 



, . • • y 



N 



A is dimensional as a 50 x 50 array. 



B - the elements of the column vectors b etc., stored in the 



fora (Bjj) 



i = 1, 2, ..., N 
j 1, 2, ..., M 



B is dimensional as a 50 x 60 array. 

X - the components of the x vectors stored in the form (Xjj) 



i - 1, 2, ..., N 
j = l, 2, ..., M 



N 
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X is dimensional as a 50 x 60 array, 

KER - error flag (See Error Return) 

2. Space Required: 

3. Temporary or Common Storage Required 

4. Error Return: Argument KER 

1, indicates no errors 

2, indicates that matrix is singular or nearly singular, 

5. Accuracy: Not applicable. 

D. METHOD 

This subroutine is an adaptation of the program F2 VTEX UNEQN. 
Please consult the write-up of that program for the mathematical method. 
The subroutine gives no printout. 



/ 
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APPENDIX D 



DESCRIPTION OF SUBROUTINE DRAW 

A. IDENTIFICATION: 

TITLE: General Graph Output Subroutine 

CO-OP ID: J7-NPS-DRAW 

CATEGORY: Output for Off-Line Plotting 

PROGRAMMER: J. R. Ward 

DATE: February 1964; REVISED JUNE 1965 

B. PURPOSE: 

This subroutine, when provided with the necessary information, generates 
a magnetic tape in the proper format for subsequent off-line graph plotting 
using the CDC 160 or 160A GRAFPLOT program and a CalComp 165 Plotter 
(see references 1 and 2). Provision is made for curve drawing and point 
plotting, automatic scaling, graph titling and axis annotating. An 
attempt was made to provide a considerable amount of flexibility, at the 
expense, necessarily, of a relatively large number of arguments and a 
rather high memory requirement. 

C. USAGE: 

1. Definitions: 

In what follows the word '‘graph 1 * will be taken to mean one piece 
or frame of graph paper on which there may be plotted one or more curves 
and/or sets of points. A "curve 1 * will mean a. continuous line generated 
by the sequence of straight lines joining successive points of the set 
defining the curve. A "point plot" will describe the representation of 
a succession of points by means of symbols (such as a cross) on the graph. 
The points are not connected in a point plot. 

2. Calling Arguments: 

All necessary information is transferred to DRAW through the calling 
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arguments. The call statement is: CALL DRAW ( NUMPTS, X,Y,MODCURV, 

I TYPE, LABEL, ITITLE, EXSCALE, YSCALE, IXUP, IYRIGHT,MODEXAX,MODEYAX, IWIDE, 
I HIGH , IGRID , LAST ) 



It is important to realize that one and only one curve or set of points 
is plotted each time DRAW is called. However, it is possible to call 
DRAW repeatedly if several curves and/or sets of points are wanted on 
one graph. 

The calling arguments are as follows: 

a, NUMPTS: The number of points defining a curve (2 ^ NUMPTS 

< 900), or the number of points to be point plotted 
(2 < NUMPTS < 30). 



b. X: 



c. I: 
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d. MODCURV: 



= 0 



= 1 



- 2 

= 3 



The array of X-ordinates ( ^X^l - 10^ for i=l, 2, . . . , 
NUMPTS). X must be dimensioned at least equal to 
NUMPTS in the calling program. 

All points will be considered to have the same 
X-ordinate if | X^y - X^ n | < 10”^. The common 
value will be put equal to zero if I^Wl ^ K>-97. 

The array of Y-ordinates, with properties corresponding 
to the .X-ordinates, above. Y must be dimensioned at 
least equal to NUMPTS in the calling program. 

Controls the number of curves, and/or sets of points 
on one graph: 

This is the only curve, or set of points, to be plotted 
on this graph. 

This is the first of two or more curves, and/or sets of 
points, to be plotted on this graph. 

This is an intermediate curve, or set of points. 

This is the last curve, or set of points, for this graph. 
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e. ITYPE 



« 0 

= 1 
- 2 
* 3 
“ 4 
“ 5 

f . LABBL: 



Controls the type of plot (i.e., curve or point plot): 
This set of points is to be represented by a curve . 
These points are to be plotted with a cross (x). 

These points are to be plotted with a plus (+). 

These points are to be plotted with a square ( □ )• 

These points are to be plotted with a diamond ( O ). 

These points are to be plotted with a triangle (A ). 

This is a Hollerith curve or point identifier. If a 
curve is being drawn, LABEL must have 4 characters 
(including any blanks), and these will be reproduced 
beside the end of the curve. This argument can be set 
in the calling program by a statement such as 
LABEL = 4H ONE ( A = blank) 

or LABEL - 4H1234 

or LABEL = 4H AAAA The latter must be used when 

no label is wanted. 

If a set of points is being plotted, LABEL is an 8- 
character identifier. The first 4 characters will be 
reproduced beside the first point, and the last four 
characters will be reproduced alongside the last point. 
This argument can be set by statements such as 
LABEL = 8HFRSTLAST 
or LABEL = 8H ONE 

or LABEL = 8H A 0NE A 123 

or LABEL *» 8H AAAAAAAA The latter must be used if 
no labels are wanted. 

The above arguments, a. through f. (and q.), have 
meaning every time DRAW is called. On the other hand. 
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g. I TITLE: 



h. EXSCALE: 



f. YSCALE: 
j. IXUP: 



k. IYRIGHT: 



1. MODEXAX: 



= 0 



the remaining arguments, g. through p., have no 
meaning except when MODCURV = 0 or 1, 

An array of twelve 8-character Hollerith words, the 
first six of which will form the first title line, and 
the last six the second. The array must be dimensioned 
12 in the calling program , must contain the user's job 
identification, and must have unwanted characters set 
to blank. For example: 

DO II = 1,12 
1 ITITLE(I) = 8H 

ITITLE(l) = 8 H a SMITH, a 
ITITLE(2) = 8HJ. A J. AAA 
ITITLE(7) = 8H a TESTIT. 

X-scale in units per inch (10”99 < EXSCALE < 1099). 
EXSCALE will always be rounded off to one figure sig- 
nificance. If EXSCALE = 0, the X-scale will be computed 
by DRAW. This is called auto-scale. 

Y-scale in units per inch, with properties corresponing 
to those of EXSCALE. 

Distance, in inches, of the X-axis from the bottom of 
the graph (0 ^ IXUP 5 IHIGH). This argument will be 
ignored unless MODEXAX = 2, see below. 

Distance, in inches, of the Y-axis from the left of the 
graph (0 — IYRIGHT — IWIDE). This will be ignored unless 
MODEYAX =2, see below. 

Determines the mode of the X-axis location: 

The X-axis will be located automatically by DRAW, with 
the origin of Y on the graph. 
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= 1 


The X-axis vdll be automatically located by DRAW, 
with the origin of Y removed (in one’s imagination) an 
integer number of inches above or below the graph, if 
this is appropriate. This option can be used only if 
the Y-scaling is automatic (YSCALE = 0). 


= 2 


The X-axis location will be as specified by IXUP. 


m. MODEYAX: 


Peteraines the mode of Y-axis. location in the same way 
as MODEXAX, above, governs the X-axis location. 


n. IWIDE: 


Width of graph in inches (l < IWIDE £ 9). If IWIDE 
is out of this range, a value of 8 will be assumed. 


o. IHIGH: 


Height of graph in inches (l < IHIGH < 15). If IHIGH 
is out of this range, a value of 8 will be assumed. 


p. IGRID: 


If IGRID = 1, a 1 M x 1" grid will be superimposed on 
the graph. This is useful only if plain paper is used 
on the CalComp Plotter. 


q. LAST: 


Indicates to the calling program whether the previous 
plot was completed successfully. The codes are: 


= 0 


Last plot was completed successfully. 


= 1 


Last plot was not completed successfully. 


* 2 


Last plot was not completed successfully, and no further 
graph output will be attempted until DRAW is next 
' entered with MODCURV = 1 or 0. 


- 3 


An attempt was made to enter DRAW with MODCURV / 1 or 
0 while the error lockout was set. 

This argument must always be a variable name and never 
a number in the call statement. 



3. Notes and Comments: 

a. The graph scales-and, if MODEXAX = 1 and/or MODEYAX = 1, the 
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amounts of origin offset - are always output as part of the graph title. 

b. Each time a graph is completed, a message to this effect is 
printed on both the standard output and the console typewriter. 

c. There are internal checks of the input to DRAW to prevent incorrect 
use. If an input error is detected, an attempt will be made, where 
possible, to complete the plot. If an argument is "corrected* 1 in this 
process, the user will be so informed on the standard output. If it is 
not possible to complete the plot, the user will be informed of the reason 
by a message on the standard output. 

d. If part or all of a curve would fall more than 0.6" laterally 
beyond the ends of the X-axis, or 0.7" vertically beyond the ends of -the 
Y-axis, the X and/or Y ordinates will be limited so that the curve will 

I typically become a line along part or all of the boundary of the graph 
as here defined. 

e. If one or more points of a point plot would fall outside the graph 
area, the plot of that point, or poirfts will be inhibited. • The number of 
such points will be reported to the user on the standard output. 

f. It should be pointed out that the X and Y scaling and axis locating 
processes are entirely independent, so that, for example, X might be 
auto-scaled, while the Y-scale is specified. At the same time the X-axis 
might be located automatically, while the Y-axis location is specified. 

g. It must be remembered that the scales and axis locations of a 
multi-plot graph are set when DRAW is called for the first time (with 
MODCURV = l). Thus the user must attempt, at that time, to achieve 
scaling and axis location which will be appropriate to all the plots he 
intends to make on the one graph. Particularly if the automatic features 
of DRAW are selected, foresight will be demanded of the user in this 
respect. 
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4. Auto-Scale Properties: 

The scale factor is chosen from amongst the values 1,2, or 5 units 

per inch, or some power of 10 times one of those values. A curve, or 

set of points which is plotted with auto-scale will normally lie entirely 

within the graph area as defined in 3.d., above. The only exception may 

occur if an axis is placed, by the user, along one edge of the graph 

(e.g., IXUP = 0, MODEXAX = 2). In such a case, points "outside** the axis 

are not considered in the selection of a scale factor (e.g., negative 

do not affect the choice of scale when IXUP = 0). If automatic axis 

location as well as auto-scale is selected, the plot, if it does not fill 

the graph area, will be placed as far as possible towards the bottom-left 

of the graph area consistant with the fact that the axes can be set in 
« 

increments of 1" only. 

5. Space Required: 3960 cells (excluding the input arrays). 

6. Temporary Storage: None. 

7. Error Print-Outs: There are a large number of possible error print- 

outs. These are all self-explanatory. 

8. Error Returns: All error returns are preceded by self-explanatory 

error printouts. An error indication is transferred back to the calling 
program through the argument LAST. 

9. Error Stops: None. 

10. Tape Mountings: Logical Tape #8 will receive the binary graph output. 

The standard monitor output will receive the messages to the user. 

11. Output Format: The format of the binary graph output records on 

magnetic tape is described in reference 1. The only difference is that 
in this program the interpolation option is by-passed (Set to zero in 
the graph output record). See reference 2. 

12. Selective Jump and Stop Settings: None. 
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13 • 1'i. fling: Variable, depending upon the number of points and the 

options chosen. Typically less than one second per curve or point plot. 

14. Accuracy: The accuracy of results is equal to the resolution of 

the CalComp Plotter, that is, 0.01** in both the X and Y directions. 

15. Equipment Configuration: CDC 1604 with FORTRAN 60 compiler and 

Library. A CDC 160 or 160A with CalComp I65 Plotter is needed for the 
off-line plotting using the appropriate GRAFPLOT program. 

D. REFERENCES: 

1. Weir, Maurice D., Spritzer, Milton and Mcllhenny, D.W., ”160-A 

Graph Plot Program,* Ident *B001, SWAP Library, 15 August 1962. 

2. Hogg, R. L. and Glover, D.C., ”160 Grafplot Routine”, Writeup 

Available from Computer Facility, U. S. Naval Postgraduate 
School, 1 April 1963. 

3. U. S. Naval Postgraduate School, Thesis, ’’Control System 

Programming, Remote Computing and Data Display,” by Robert Lee 
Hogg and Dennis C. Glover, 1963. 



N.B, References 1 and 2 are included in Reference 3 as 
Enclosures 2 and 1, respectively. 
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APPENDIX E 



PLOTS OF DOCK TRIAL DATA 

Scale of stress has been omitted for reasons of security 
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FIGURE E-l 

Polar plot of On for rosette i^Lng #1 
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FIGURE E-2 

Polar plot of 'tx.y for rosette ring #1 
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"FIGURE E-3 



Polar plot of 0 *kj for rosette ring #2 

86 



1.0 

CS> 




CO 

< 5 > 



or, 

< 5 > 



FIGURE E-4 

Polar plot of for rosette ring # 2 
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FIGURE E-5 

Polar plot of On for rosette ring #3 
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FIGURE E-6 

y JWter plot pf 'key for rosette ring #3" 
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FIGURE E-7 

Polar plot of (Tm for rosette ring #4 
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FIGURE E-9 

•.« Polar plot of CTn for rosette ring #5 
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FIGURE E-10 

Polar .plot of 't*Y for rosette ring #5 
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FIGURE E-12 

Polar plot of 'txy for rosette ring #6 
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FIGURE E-13 

Polar plot of CT n for rosette ring #7 
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FIGURE E-I4 

Polar plot of 'txy for rosette ring #7 
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FIGURE E-15 



Polar plot of On for rosette ring #8 
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FIGURE E-16 



Polar plot of *txy for rosette ring #8 
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